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In a recent paper, Teo and Kane proposed a 3D model in which the defects support Majorana fermion zero 
modes. They argued that exchanging and twisting these defects would implement a set TZ of unitary trans- 
formations on the zero mode Hilbert space which is a 'ghostly' recollection of the action of the braid group 
on Ising anyons in 2D. In this paper, we find the group 72n which governs the statistics of these defects by 
analyzing the topology of the space K2n of configurations of 2n defects in a slowly spatially- varying gapped 
free fermion Hamiltonian: Tin = 7ri(^2n). We find that the group Tin = Z x T-fn, where the 'ribbon per- 
mutation group' Tin is a mild enhancement of the permutation group Sm'- Tin = Z2 x E(('Li)'^'^ x: Sin)- 
Here, Eiili)^^ xi Si„) is the 'even part' of {"LiY" x Sin, namely those elements for which the total parity 
of the element in iJjif" added to the parity of the permutation is even. Surprisingly, TZ is only a projective 
representation of Tin, a possibility proposed by Wilczek. Thus, Teo and Kane's defects realize 'Projective Rib- 
bon Permutation Statistics', which we show to be consistent with locality. We extend this phenomenon to other 
dimensions, co-dimensions, and symmetry classes. Since it is an essential input for our calculation, we review 
the topological classification of gapped free fermion systems and its relation to Bott periodicity. 
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PACS numbers: 

I. INTRODUCTION 

In two dimensions, the configuration space of n point- 
like particles C^^ is multiply-connected. Its first homotopy 
group, or fundamental group, is the n-particle braid group, 
7ri(C^^) = Bn- The braid group Bn is generated by counter- 
clockwise exchanges ai of the i* and (i + 1)* particles satis- 
fying the defining relations: 



cr,:CTi 



(T,(7,: 



for |«- jl > 2 



CiCTj+lCi — fTi+lCTi (Ti+l 



for 1 < z < n - 2 (1) 



This is an infinite group, even for only two particles, since 
((Ti)™ is a non-trivial element of the group for any tti > 0. In 
fact, even if we consider distinguishable particles, the result- 
ing group, called the 'pure Braid group' is non-trivial. (For 
two particles, the pure braid group consists of all even powers 

OfCTl.) 

In quantum mechanics, the equation 7ri(C^^) = Bn opens 
the door to the possibility of anyons''^. Higher-dimensional 
representations of the braid group give rise to non-Abelian 
anyons^"^. There has recently been intense effort directed 
towards observing non-Abelian anyons due, in part, to their 
potential use for fault-tolerant quantum computation*"'^. One 
of the simplest models of non-Abelian anyons is called Ising 
anyons. They arise in theoretical models of the v — 5/2 frac- 
tional quantum Hall state*^"' ' (see also Ref. 12), chiral p-wave 
superconductors'^''', a solvable model of spins on the hon- 
eycomb lattice'^, and interfaces between superconductors and 
either 3D topological insulators '*" or spin-polarized semicon- 
ductors with strong spin-obrit coupling'^. A special feature 
of Ising anyons, which makes them relatively simple and con- 
nects them to BCS superconductivity, is that they can be un- 
derstood in a free fermion picture. 

A collection of 2n Ising anyons has a 2"^'^ -dimensional 



Hilbert space (assuming fixed boundary condition). This can 
be understood in terms of 2n Majorana fermion operators 

7j = 7J, i = 1, 2, . . . , n, one associated to each Ising anyon, 
satisfying the anticommutation rules 



{l^^j} = 2(5,^ 



(2) 



The Hilbert space of 2n Ising anyons with fixed boundary con- 
dition furnishes a representation of this Clifford algebra; by 
restricting to fixed boundary condition, we obtain a represen- 
tation of products of an even number of 7 matrices, which has 
minimal dimension 2""'^. When the i* and {i + 1)* anyons 
are exchanged in a counter-clockwise manner, a state of the 
system is transformed according to the action of 



P{<^t 



gi-n'/S g-'r7i7i + i/4 



(3) 



(There is a variant of Ising anyons, associated with SU(2)2 
Chern-Simons theory, for which the phase factor e^^/^ is re- 
placed by e~"^/*. In the fractional quantum Hall effect, Ising 
anyons are tensored with Abelian anyons to form more com- 
plicated models with more particle species; the phase factor 
depends on the model.) A key property, essential for appli- 
cations to quantum computing, is that a pair of Ising anyons 
forms a two-state system. The two states correspond to the 
two eigenvalues ±1 of 7^7,. No local degree of freedom can 
be associated with each anyon; if we insisted on doing so, we 
would have to say that each Ising anyon has a/2 internal states. 
In superconducting contexts, the 7jS are the Bogoliubov-de 
Gennes operators for zero-energy modes (or, simply, 'zero 
modes') in vortex cores; the vortices themselves are Ising 
anyons if they possess a single such zero mode 7^. Although 
the Hilbert space is non-local in the sense that it cannot be de- 
composed into the tensor product of local Hilbert spaces as- 
sociated with each anyon, the system is perfectly compatible 



with locality and arises in local lattice models and quantum 
field theories. 

In three or more dimensions, the configuration space of 
n point-like particles is simply-connected if the particles are 
distinguishable. If the particles are indistinguishable, it is 
multiply-connected, 7ri(C^^) = 5„. The generators of the 
permutation group satisfy the relations (1) and one more, 
erf = 1. As a result of this last relation, the permutation group 
is finite. The one-dimensional representations of Sn corre- 
spond to bosons and fermions. One might have hoped that 
higher-dimensional representations of Sn would give rise to 
interesting 3D analogues of non-Abelian anyons. However, 
this is not the case, as shown in Ref. 18,19: any higher- 
dimensional representation of 5„ which is compatible with 
locality can be decomposed into the tensor product of local 
Hilbert spaces associated with each particle. For instance, 
suppose we had 2n spin- 1/2 particles but ignored their spin 
values. Then we would have 2^" states which would trans- 
form into each other under permutations. Clearly, if we dis- 
covered such a system, we would simply conclude that we 
were missing some quantum number and set about trying to 
measure it. This would simply lead us back to bosons and 
fermions with additional quantum numbers. (The color quan- 
tum number of quarks was conjectured by essentially this kind 
of reasoning.) The quantum information contained in these 
2^" states would not have any special protection. 

The preceding considerations point to the following ten- 
sion. The Clifford algebra (2) of Majorana fermion zero 
modes is not special to two dimensions. One could imagine 
a three (or higher) dimensional system with topological de- 
fects supporting such zero modes. But the Hilbert space of 
these topological defects would be 2"^^-dimensional, which 
manifestly cannot be decomposed into the tensor product 
of local Hilbert spaces associated with each particle, seem- 
ingly in contradiction with the results of Refs. 18,19 on 
higher-dimensional representations of the permutation group 
described above. However, as long as no one had a three 
or higher dimensional system in hand with topological de- 
fects supporting Majorana fermion zero modes, one could, 
perhaps, sweep this worry under the rug. Recently, however, 
Teo and Kane^" have shown that a 3D system which is simul- 
taneously a superconductor and a topological insulator^ '"^"^ 
(which, in many but not all examples, is arranged by forming 
superconductor-topological insulator heterostructures) sup- 
ports Majorana zero modes at point-like topological defects. 

To make matters worse, Teo and Kane^*' further showed 
that exchanging these defects enacts unitary operations on this 
2"- 1 -dimensional Hilbert space which are essentially equal to 
(3). But we know that these unitary matrices form a represen- 
tation of the braid group, which is not the relevant group in 
3D. One would naively expect that the relevant group is the 
permutation group, but S'„ has no such representation (and 
even if it did, its use in this context would contradict locality, 
according to Ref. 18,19 and arguments in Ref. 25). So this 
begs the question: what is the group 72n for which Teo and 
Kane's unitary transformations form a representation? 

With the answer to this question in hand, we could address 
questions such as the following. We know that a 3D incarna- 



tion of Ising anyons is one possible representation of 72ri; is a 
3D version of other anyons another representation of Tin? 

Attempts to generalize the braiding of anyons to higher di- 
mensions sometimes start with extended objects, whose con- 
figuration space may have fundamental group which is richer 
than the permutation group. Obviously, if one has line-like de- 
fects in 3D which are all oriented in the same direction, then 
one is essentially back to the 2D situation governed by the 
braid group. This is too trivial, but it is not clear what kind 
of extended objects in higher dimensions would be the best 
starting point. What is clear, however, is that Teo and Kane's 
topological defects must really be some sort of extended ob- 
jects. This is clear from the above-noted contradiction with 
the permutation group. It also follows from the 'order param- 
eter' fields which must deform as the defects are moved, as 
we will discuss. 

In this paper, we show that Teo and Kane's defects are prop- 
erly viewed as point-like defects connected pair-wise by rib- 
bons. We call the resulting 2n-particle configuration space 
K2n- We compute its fundamental group TTi{K2n), which 
we denote by 72n and find that 72ri = Z x T-J^- Here, 
T2n is the 'ribbon permutation group', defined by T^i = 
Z2 X S((Z2)2« X 52„). The group ^((Za)^" x S2n) is a 
non-split extension of the permutation group S'2„ by Zg"^^ 
which is defined as follows: it is the subgroup of (Z2)^" x 5*2™ 
composed of those elements for which the total parity of the 
element in (Z2)^" added to the parity of the permutation is 
even. The 'ribbon permutation group' for 2n particles, by TJ^, 
is the fundamental group of the reduced space of 2n-particle 
configurations. 

Our analysis relies on the topological classification of 
gapped free fermion Hamiltonians^^-^^ - band insulators and 
superconductors - which is the setting in which Teo and 
Kane's 3D defects and their motions are defined. The start- 
ing point for this classification is reducing the problem from 
classifying gapped Hamiltonians defined on a lattice to clas- 
sifying Dirac equations with a spatially varying mass term. 
One can motivate the reduction to a Dirac equation as Teo 
and Kane do: they start from a lattice Hamiltonian and as- 
sume that the parameters in the Hamiltonian vary smoothly 
in space, so that the Hamiltonian can be described as a func- 
tion of both the momentum k and the position r. Near the 
minimum of the band gap, the Hamiltonian can be expanded 
in a Dirac equation, with a position-dependent mass term. In 
fact, Kitaev^^ has shown that the reduction to the Dirac equa- 
tion with a spatially varying mass term can be derived much 
more generally: gapped lattice Hamiltonians, even if the pa- 
rameters in the Hamiltonian do not vary smoothly in space, 
are stably equivalent to Dirac Hamiltonians with a spatially 
varying mass term. Here, equivalence of two Hamiltonians 
means that one can be smoothly deformed into the other while 
preserving locality of interactions and the spectral gap, while 
stable equivalence means that one can add additional "triv- 
ial" degrees of freedom (additional sites which have vanishing 
hopping matrix elements) to the original lattice Hamiltonian 
to obtain a Hamiltonian which is equivalent to a lattice dis- 
cretization of the Dirac Hamiltonian. 

Since this classification of Dirac Hamiltonians is essential 



for the definition of K2n, we give a self-contained review, 
following Kitaev's analysis^^. Our exposition parallels the 
discussion of Bott periodicity in Milnor's book^^. The basic 
idea is that each additional discrete symmetry which squares 
to —1 which we impose on the system is encapsulated by 
an anti- symmetric matrix which defines a complex structure 
on R^, where N/2 is the number of bands (or, equivalently, 
A^ is the number of bands of Majorana fermions). For any 
given system, these are chosen and fixed. This leads to a pro- 
gression of symmetric spaces 0(iV) -^ 0{N)/U{N/2) -^ 
U{N/2)/Sp{N/4) -^ ... as the number of such symmetries 
is increased. Following Kitaev^^, we view the Hamiltonian 
as a final anti-symmetric matrix which must be chosen (and, 
thus, put almost on the same footing as the symmetries); it is 
defined by a choice of an arbitrary point in the next symmet- 
ric space in the progression. The space of such Hamiltonians 
is topologically-equivalent to that symmetric space. However, 
as the spatial dimension is increased, 7-matrices squaring to 
+1 must be chosen in order to expand the Hamiltonian in the 
form of the Dirac equation in the vicinity of a minimum of 
the band gap. These halve the dimension of subspaces of M^ 
by separating it into their +1 and —1 eigenspaces and thereby 
lead to the opposite progression of symmetric spaces. Thus, 
taking into account both the symmetries of the system and 
the spatial dimension, we conclude that the space of gapped 
Hamiltonians with no symmetries in d = 3 is topologically 
equivalent to {J{N)/0{N). (However, by the preceding con- 
siderations, the same symmetric space also, for instance, clas- 
sifies systems with time-reversal symmetry in d = 4.) All 
such Hamiltonians can be continuously deformed into each 
other without closing the gap, 7ro(U(7V)/0(iV)) = 0. How- 
ever, there are topologically-stable point-like defects classi- 
fied by 7r2(U(iV)/0(iV)) = Z2. These are the defects whose 
multi-defect configuration space we study in order to see what 
happens when they are exchanged. 

The second key ingredient in our analysis is 1950's-vintage 
homotopy theory, which we use to compute TTi{K2n)- We 
apply the Pontryagin-Thom construction to show that K2n, 
which includes not only the particle locations but also the full 
field configuration around the particles (i.e. the way in which 
the gapped free fermion Hamiltonian of the system explores 
U{N)/0{N)), is topologically-equivalent to a much simpler 
space, namely point-like defects connected pair-wise by rib- 
bons. In order to then calculate Tri(K2n), we rely on the long 
exact sequence of homotopy groups 



7r,{E) ^ 7r,(B) ^ 7r,_i(F) ^ 7r,_i(S) 



(4) 



associated to a fibration defined by F ^ i? — > S. (In an 
exact sequence, the kernel of each map is equal to the image 
of the previous map.) This exact sequence may be familiar to 
some readers from Mermin's review of the topological theory 
of defects^^, where a symmetry associated with the group G 
is spontaneously broken to H, thereby leading to topological 
defects classified by homotopy groups 7r„ (G/iJ). These can 
be computed by (4) with E = G,F = H,B = G/H, e.g. if 
7ri(G) = 7ro(G) = 0, then 7ri(G/iJ) = na{H). 

The ribbon permutation group is a rather weak enhance- 



ment of the permutation group and, indeed, we conclude that 
Teo and Kane's unitary operations are not a representation of 
the ribbon permutation group. However, they are a projective 
representation of the ribbon permutation group. In a projective 
representation, the group multiplication rule is only respected 
up to a phase, a possibility allowed in quantum mechanics. A 
representation p (sometimes called a linear representation) of 
some group G is a map from the group to the group of lin- 
ear transformations of some vector space such that the group 
multiplication law is reproduced: 



p{gh) = pig) ■ p{h) 



(5) 



if .g, h E G. Particle statistics arising as a projective represen- 
tation of some group realizes a proposal of Wilczek's-'", albeit 
for the ribbon permutation group rather than the permutation 
group itself This difference allows us to sidestep a criticism 
of Read^^ based on locality, which Teo and Kane's projective 
representation respects. The group (Z2)^"^^ is generated by 
2n — 1 generators xi,X2, ■■■, a;2„-i satisfying 



ijj l"^ t 



1 



(6) 



However, the projective representation of (Z2)^"^^, which 
gives a subgroup of Teo and Kane's transformations, is an or- 
dinary linear representation of a Z2-central extension, called 
the extra special group i?2ri-i- 






= 1 



XjXi for \i — j\>2 

1 



(7) 



Here, z generates the central extension, which we may take to 
be z — —1. The operations generated by the x^s were dubbed 
'braidless operations' by Teo and Kane^" because they could 
be enacted without moving the defects. While these opera- 
tions form an Abelian subgroup of 72n, their representation 
on the Majorana zero mode Hilbert space is not Abelian - two 
such operations which twist the same defect anti-commute 
(e.g. Xi and x^+i). 

The remaining sections of this paper will be as follows. In 
Section II, we rederive Teo and Kane's zero modes and unitary 
transformations by simple pictorial and counting arguments 
in a 'strong-coupling' limit of their model. In Section III, we 
review the topological classification of free-fermion Hamilto- 
nians, including topological insulators and superconductors. 
From this classification, we obtain the classifying space rel- 
evant to Teo and Kane's model and, in turn, the topological 
classification of defects and their configuration space. In Sec- 
tion IV, we use a toy model to motivate a simple picture for the 
defects used by Teo and Kane and give a heuristic construction 
of the ribbon permutation group. In Section V, we give a full 
homotopy theory calculation. In Section VI, we compare the 
ribbon permutation group to Teo and Kane's unitary transfor- 
mations and conclude that the latter form a projective, rather 
than a linear, representation of the former. Finally, in Section 
VII, we review and discuss our results. Several appendices 



contain technical details. 



II. STRONG-COUPLING LIMIT OF THE TEO-KANE 
MODEL 

In this section, we present a lattice model in d dimen- 
sions which has, as its continuum limit in d = 3, the model 
discussed by Teo and Kane^". In the limit that the mass 
terms in this model are large (which can be viewed as a 
'strong-coupling' limit), a simple picture of topological de- 
fects ('hedgehogs') emerges. We show by a counting ar- 
gument that hedgehogs possess Majorana zero modes which 
evolve as the hedgehogs are adiabatically moved. This adia- 
batic evolution is the 3D non-Abelian statistics which it is the 
main purpose of this paper to explain. 

The strong coupling limit which we describe is the sim- 
plest way to derive the existence of Majorana zero modes and 
the unitary transformations of their Hilbert space which re- 
sults from exchanging them. This section does not require the 
reader to be au courant with the topological classification of 
insulators and superconductors'*"'^^. (In the next section, we 
will review that classification in order to make our exposition 
self-contained.) 

We use a hypercubic lattice in d-dimensions, with a single 
Majorana degree of freedom at each site. That is, for d ~ \, 
we use a chain, in d = 2 we use a square lattice, in rf = 3 
we use a cubic lattice, and so on. We first construct a lattice 
model whose continuum limit is the Dirac equation with 2^^- 
dimensional 7-matrices to reproduce the Dirac equation con- 
sidered by Teo and Kane; we then show how to perturb this 
model to open a mass gap. We begin by considering only 
nearest neighbor couplings. The Hamiltonian H is an anti- 
symmetric Hermitian matrix. In d = 1, we can take the linear 
chain to give a lattice model with the Dirac equation as its 



a) 



b) 



continuum limit. That is, i?j.,+i — i and H 



2j+ij 



To 



describe this state in pictures, we draw these bonds as oriented 
lines, as shown in Fig. (la), with the orientation indicating the 
sign of the bond. The continuum limit of this Hamiltonian is 
described by a Dirac equation with 2-dimensional 7 matrices. 
While this system can be described by a unit cell of a sin- 
gle site, we instead choose to describe it by a unit cell of two 
sites for convenience when adding mass terms later. In rf = 2, 
we can take a 7r-flux state to obtain the Dirac equation in the 
continuum limit. A convenient gauge to take to describe the 
TT-flux state is shown in Fig. (lb), with all the vertical bonds 
having the same orientation, and the orientation of the hor- 
izontal bonds alternating from row to row. The continuum 
limit here has 4-dimensional 7 matrices and we use a 4-site 
unit cell. 

In general, in d dimensions, we can obtain a Dirac equation 
with 2''-dimensional 7 matrices by the following iterative pro- 
cedure. Let the "vertical" direction refer to the direction of the 
d-th basis vector. Having constructed the lattice Hamiltonian 
ind—l dimensions, we stack these Hamiltonians vertically on 
top of each other, with alternating signs in each layer. Then, 
we take all the vertical bonds to be oriented in the same di- 
rection. This Hamiltonian is invariant under translation in the 
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FIG. 1 : (a) A lattice model giving the Dirac equation in d = 1. (b) A 
lattice model in d = 2. 



FIG. 2: Dimerization in d = 1. 



vertical direction by distance 2. Thus, if Hd-i is the Hamil- 
tonian in d — 1 dimensions, the Hamiltonian Hd is given by 



H, 



Hd-i 2sin(/c/2)A 
2sin(fc/2)/ -Hd-i ) 



(8) 



where / is the identity matrix and k is the momentum in the 
vertical direction. Near fc = 0, this is 



Hd ~ H, 



d-l 



(9) 



This iterative construction corresponds to an iterative con- 
struction of 7-matrices. Having constructed d — 1 different 
2''^ 1 -dimensional 7-matrices 7^ , . . . , 7^_ i , we construct d dif- 
ferent 2''-dimensional 7-matrices, 7]^, ..., 7^, by 7^ = 7^ ® (Jz 
for z = 1, ..., d — 1, and j^ = I (E) (Tx- 

In one dimension, dimerization of bonds corresponds to al- 
ternately strengthening and weakening the bonds as shown in 
Fig. (2). In two dimensions, we can dimerize in either the 
horizontal or vertical directions. In d-dimensions, we have d 
different directions to dimerize. Dimerizing in the "vertical" 
direction gives, instead of (9), the result 



Hd 



Hd-i 



rud ®(Ty, 



(10) 



where md is the dimerization strength. This corresponds to 
an iterative construction of mass matrices. Mi, as follows. In 
one dimension, we have Mi = idy. Given d — 1 different 
mass matrices in d — 1 dimensions. Mi, we construct Mi in 
d-dimensions by Mi = Mi (E) a^, for i = l...d — 1, and 
Md ^iI(E)(Ty. 

If the dimerization is non-zero, and constant, we can in- 
crease the dimerization strength without closing the gap un- 
til a strong coupling limit is reached. In one dimension, by 
increasing the dimerization strength, we eventually reach a 
fully dimerized configuration, in which each site has one non- 
vanishing bond connected to it. In two or more dimensions, 
the dimerization can be a combination of dimerization in dif- 
ferent directions. However, if the dimerization is completely 
in one direction, for example the vertical direction, we in- 
crease the dimerization strength until the vertical bonds are 
fully dimerized. Simultaneously, we reduce the strength of 
the other bonds to zero without closing the gap. This is again 



a fully dimerized state, the columnar state, with each site hav- 
ing one non-vanishing bond. Any configuration with uniform, 
small dimerization can be deformed into this pattern with- 
out closing the gap by rotating the direction of dimerization, 
increasing the strength of dimerization, and then setting the 
bonds in the other directions to zero. 

It is important to understand that the ability to reach such 
a strong coupling limits depends on the perturbation of the 
Dirac equation that we consider; for dimerization, it is possi- 
ble to reach a strong coupling limit, while if we had instead 
chosen to open a mass gap by adding, for example, diago- 
nal bonds with imaginary coupling to the two-dimensional 
Dirac equations, we would open a mass gap by perturbing 
the Hamiltonian with the term 17172, and such a perturba- 
tion cannot be continued to the strong coupling limit due to 
topological obstruction. 

Further, if the dimerization is non-uniform then it may not 
be possible to reach a fully dimerized state without having 
defect sites. Consider the configurations in Fig. (3a) in d = 1 
and in Fig. (3b) in d = 2. These are the strong coupling limits 
of the hedgehog configuration, and each contains a zero mode, 
a single unpaired site. This is one of the central results of the 
strong-coupling limit: topological defects have unpaired sites 
which, in turn, support Majorana zero modes. 

Such strong-coupling hedgehog configurations can be con- 
structed by the following iterative process in any dimension d. 
Let x^ correspond to the coordinate in the vertical direction. 
For Xd > 0, stack d — 1 -dimensional hedgehog configura- 
tions. Along the half-line given by a:;^ > and x^ = for 
1 < i < d — 1, arrange vertical bonds, oriented to connect 
the site with x^ = 2fc — 1 to that with x^ — 2k, for fc > 1. 
Along the lower half plane, given by Xd < 0, arrange verti- 
cal bonds oriented to connect a site with Xd — —{2k — 1) to 
that with —2k, for k > 1. This procedure gives the d — 2 
hedgehog in Fig. (3b) from the d = 1 hedgehog in Fig. (3a), 
and gives a strong coupling limit of the Teo-Kane hedgehog 
in d = 3. That is, the Teo-Kane hedgehog can be deformed 
into this configuration, without closing the gap. 

So long as we consider only nearest-neighbor bonds, there 
is an integer index v describing different dimerization pat- 
terns in the strong-coupling limit. This index, which is present 
in any dimension, arises from the sublattice symmetry of the 
system, and is closely-related to the U(l) symmetry of dimer 
models of spin systems^'. Label the two sublattices by A and 
B. Consider any set of sites, such that every site in that set has 
exactly one bond connected to it. (Recall that, in the strong 
coupling limit, every bond has strength or 1 and every site 
has exactly one bond connected to it, except for defect sites.) 
Then, the number of bonds going from A sites in this set to 
B sites outside the set is exactly equal to the number of bonds 
going from B sites in this set to A sites outside the set. On 
the other hand, if there are defect sites in the set, then this rule 
is broken. Consider the region defined by the dashed line in 
Fig. (4a). We define the "flux" crossing the dashed line to be 
the number of bonds crossing that boundary which leave start- 
ing on an A site, minus the number which leave starting on a 
B site. The flux is non-zero in this case, but is unchanging as 
we increase the size of the region. This flux is the index jy. By 
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b) 
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FIG. 3: (a) A one-dimensional hedgehog, (b) A two-dimensional 
hedgehog. 



the argument given above for the existence of zero modes, i^ 
computed for any region is equal to the number of Majorana 
zero modes contained within the region. 

The index v can be defined beyond the strong-coupling 
limit. Consider, for the sake of concreteness, d = 3. There 
are 3 possible dimerizations, one for each dimension, as we 
concluded in Eq. 10. In weak-coupling, the square of the gap 
is equal to the sum of the squares of the dimerizations. Thus, 
if we assume a fixed gap, we can model these dimerizations 
by a unit vector. The integer index discussed above is simply 
the total winding number of this unit vector on the boundary 
of any region. 

However, once diagonal bonds are allowed, the integer in- 
dex ly no longer counts zero modes. Instead, there is a Z2 
index, equal to z/(mod2) which counts zero modes modulo 2. 
To see this in the strong-coupling limit, consider the configu- 
ration in Fig. (4b). This is a configuration with i^ = 2 but no 
Majorana zero modes. However, a i^ = 1 configuration must 
still have a zero mode and, thus, any configuration with odd i^ 
must have at least one zero mode. 

In Fig. (4), we have chosen to orient the bonds from A to B 
sublattice to make it easier to compute z/. However, the ly and 
its residue modulo 2, defined above are independent of the ori- 
entation of the bonds (which indicate the sign of terms in the 
Hamiltonian) and depend only on which sites are connected 
by bonds (which indicate which terms in the Hamiltonian are 
non-vanishing). 

The i^(mod2) with diagonal bonds is the same as Kitaev's 
"Majorana number"'^. We can use this to show the exis- 
tence of zero modes in the Teo-Kane hedgehog even outside 
the strong-coupling limit. Consider a hedgehog configuration. 
Outside some large distance R from the center of the hedge- 
hog, deform to the strong coupling limit without closing the 
gap. Then, outside a distance R, we can count i^(mod2) by 
counting bonds leaving the region and we find a nonvanishing 
result relative to a reference configuration: if there are an even 
number of sites in the region then there are an odd number of 



b) 



.• 



I I 



I I 



FIG. 4: (a) Defect acting as source of (7(l)-flux. Bonds are oriented 
from A to B sublattice. There is a net flux of one leaving the region 
defined by the dashed line, (b) Configuration with diagonal bond 
added, indicated by the undirected line connecting the two circles; 
either orientation of this line, corresponding to different choices of 
the sign of the term in the Hamiltonian, would lead to the same result. 
There is a net of flux of two leaving the region defined by the dashed 
line. 



bonds leaving in a hedgehog configuration, and if there are an 
odd number of sites then there are an even number of bonds 
leaving. However, since this implies a nonvanishing Majorana 
number, there must be a zero mode inside the region, regard- 
less of what the Hamiltonian inside is. We note that this is a 
highly non-trivial result in the weak-coupling limit, where the 
addition of weak diagonal bonds, all oriented the same direc- 
tion, to the configuration of Fig. (lb) corresponds to adding 
the term ijiJ2 to the Hamiltonian in rf = 2. By the argument 
given above, even this Hamiltonian has a zero mode in the 
presence of a defect with non-zero z/(mod2). 

Given any two zero modes, corresponding to defect sites in 
the strong coupling limit, we can identify a string of sites con- 
necting them. If we have a pair of defect sites on opposite 
sublattices, corresponding to opposite hedgehogs, then one 
particular string corresponds to the north pole of the order 
parameter, as in Fig. (5a). However, we can simply choose 
any arbitrary string. Let 7^,7, be the Majorana operators at 
the two defect sites. The operation 7^ ^- —7^,7, — > —7, 
can be implemented as follows. We begin with an adiabatic 
operation on one of the defect sites and the nearest 2 sites 
on the line. The Hamiltonian on those three sites is an anti- 
symmetric, Hermitian matrix. That is, it corresponds to an 
oriented plane in three dimensions. We can adiabatically per- 
form orthogonal rotations of this plane. Thus, by rotating by 
TT in the plane corresponding to the defect site and the first 
site on the string, we can change the sign of the mode on the 
defect and the orientation of the bond, as shown in Fig. (5b). 
This rotation is an adiabatic transformation of the three site 
Hamiltonian 



ism{e) 

i cos(0) 

-isin{9) —icos{9) 



(11) 



along the path = ^ tt. We then perform rotations on con- 
secutive triples of sites along the defect line, which changes 
the orientation of pairs of neighboring bonds, arriving at the 
configuration in Fig. (5c). Finally, we rotate by n in the plane 
containing the other defect site and the last site. This returns 



FIG. 5: (a) Pairs of defects connected by a string, (b) First rotation 
applied to the configuration in (a) Open circle replaces filled circle 
to indicate sign change of the Majorana mode on the site. (c)After 
rotating along the string, (d) Rotating the last site and restoring the 
string to its original configuration 



the system to the original configuration, having effected the 
desired operation. 

Since we only consider adiabatic transformation, we can 
only perform orthogonal rotations with unit determinant. 
Thus, any transformation which swaps two defects and returns 
the bonds to their original configuration, must change the sign 
of one of the zero modes: 7^ -^ 7 , 7 — > — 7^. Indeed, any or- 
thogonal transformation with determinant equal to minus one 
would change the sign of the fermion parity in the system, as 
the fermion parity operator is equal to the product of the ji 
operators. 

We used the ability to change the orientation of a pair of 
bonds in this construction. The fact that one can only change 
the orientation of bonds in pairs, and not the orientation of a 
single bond, is related to a global Z2 invariant: the Hamilto- 
nian is an anti-symmetric matrix and the sign of its Pfaffian 
cannot be changed without closing the gap. Changing the di- 
rection of a single bond changes the sign of this Pfaffian and 
so is not possible. 

The above discussion left open the question of which zero 
changes its sign, i.e. is the effect of the exchange 7^ -^ 
lj-,lj -^ -li or 7j -^ ~7j'7i ~^ li^ The answer is 
that it depends on how the bonds are returned to their orig- 
inal configuration after the exchange is completed (which is 
a clue that the defects must be understood as extended ob- 
jects, not point-like ones). For the bonds to be restored, one 
of the defects must be rotated by 27r; the corresponding zero 
mode acquires a minus sign. We will discuss this in greater 
detail in a later section. The salient point here is that the 
effect of an exchange is a unitary transformation generated 
by the operator e '''•^j'^. This is reminiscent of the repre- 
sentation of braid group generators for non-Abelian quasi- 
particles in the quantum Hall effect'' and vortices in chiral 
p-wave superconductors '"*, namely the braid group represen- 
tation realized by Ising anyons^. But, of course, in 3D the 
braid group is not relevant, and the permutation group, which 
is associated with point-like particles in d > 2, does not have 
non-trivial higher-dimensional representations consistent with 
locality'**''^. As noted in the introduction, this begs the ques- 



tion: what group are the unitary matrices e 
ing? 



±7r7.7./4 



represent- 



III. TOPOLOGICAL CLASSIFICATION OF GAPPED 
FREE FERMION HAMILTONIANS 

A. Setup of the Problem 

In this section, we will briefly review the topological classi- 
fication of translationally-invariant or slowly spatially-varying 
free-fermion Hamiltonians following Kitaev's analysis in Ref. 
27. (For a different perspective, see Schnyder et al. 's approach 
in Ref. 26.) The 3D Hamiltonian of the previous section is a 
specific example which fits within the general scheme and, by 
implication, the 3D non-Abelian statistics which we derived at 
the end of the previous section also holds for an entire class of 
models into which it can be deformed without closing the gap. 
Our discussion will follow the logic of Milnor's treatment of 
Bott periodicity in Ref. 28. 

Consider a system of N flavors of electrons Cj (k) in d di- 
mensions. The flavor index j accounts for spin as well as the 
possibility of multiple bands. Since we will not be assuming 
charge conservation, it is convenient to express the complex 
fermion operators Cj (k) in terms of real fermionic operators 
(Majorana fermions), Cj(k) = (02^-1 (k) + ia2j(k))/2 (the 
index j now runs from 1 to 2A''). The momentum k takes 
values in the Brillouin zone, which has the topology of the d- 
dimensional torus T''. The Hamiltonian may be written in the 
form 



H 






(12) 



where, by Fermi statistics, Aij{p) = — Aji(— p). Let us 
suppose that the Hamiltonian (12) has an energy gap 2 A, by 
which we mean that its eigenvalues E^ (p) (a is an index la- 
beling the eigenvalues of H) satisfy |£^Q,(p)| > A. The basic 
question which we address in this section is the following. 
What topological obstructions prevent us from continuously 
deforming one such gapped Hamiltonian into another? 

Such an analysis can apply, as we will see, not only to free 
fermion Hamiltonians, but also to those interacting fermion 
Hamiltonians which, deep within ordered phases, are well- 
approximated by free-fermion Hamiltonians. (This can in- 
clude rather non-trivial phases such as Ising anyons, but not 
Fibonacci anyons.) In such settings, the fermions may be 
emergent fermionic quasiparticles; if the interactions between 
these quasiparticles are irrelevant in the renormalization- 
group sense, then an analysis of free-fermion Hamiltonians 
can shed light on the phase diagrams of such systems. Thus, 
the analysis of free fermion Hamiltonians is equivalent to 
the analysis of interacting fermion ground states whose low- 
energy quasiparticle excitations are free fermions. 

Let us begin by considering a few concrete examples, in 
order of increasing complexity. 



B. Zero-Dimensional Systems 

First, we analyze a zero-dimensional system which we will 
not assume to have any special symmetry. The Hamiltonian 



(12) takes the simpler form: 



(13) 



where Aij is a 2N x 2N antisymmetric matrix, Aij = —Aji. 
Any real antisymmetric matrix can be written in the form 
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where O is an orthogonal matrix and the A^s are positive. The 
eigenvalues of A come in pairs izA^; thus, the absence of 
charge conservation can also be viewed as the presence of a 
particle-hole symmetry. By assumption, A^ > A for all i. 
Clearly, we can continuously deform A^ without closing the 
gap so that A^ = A for all i. (This is usually called 'spectrum 
flattening'.) Then, we can write: 



where 



yl = A • d^JO 
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The possible choices of A^ correspond to the possible choices 
of O e 0(2Af), modulo O which commute with the matrix 
J. But the set of O e 0(2iV) satisfying O^ JO = J 
is U(A^)cO(2iV). Thus, the space of all possible zero- 
dimensional free fermionic Hamiltonians with N single- 
particle energy levels is topologically-equivalent to the sym- 
metric space 0(2A^)/U(7V). 

This can be restated in more geometrical terms as follows. 
Let us here and henceforth take units in which A = L Then 
the eigenvalues of A are ±i. If we view the 2iV x 2iV ma- 
trix A as a linear transformation on M^^, then it defines a 
complex structure. Consequently, we can view M^^ as C^ 

.^^ by a complex scalar can 
+ hAv. The set of complex 



since multiplication of w G H 
be defined as (a + i\))v = av 



structures on M is given by performing an arbitrary 0(2A^) 
rotation on a fixed complex structure, modulo the rotations of 
C^ which respect the complex structure, namely U(A^). Thus, 
once again, we conclude that the desired space of Hamiltoni- 
ans is topologically-equivalent to 0(2N)/U(N). 

What are the consequences of this equivalence? Con- 
sider the simplest case, N = 1. Then, the space of 
zero-dimensional Hamiltonians is topologically-equivalent to 
0(2)/U(l)= Z2: there are two classes of Hamiltonians, those 
in which the single fermionic level is unoccupied in the 
ground state, c^c = (1 + iaia2)/2 = 0, and those in which 
it is occupied. For larger N, 0{2N)flJ(N) is a more com- 



plicated space, but it still has two connected components, 
TTo{0{2N) /U{N)) — Z2, so that there are two classes of free 
fermion Hamiltonians, corresponding to even or odd numbers 
of occupied fermionic levels in the ground state. 

Suppose, now, that we restrict ourselves to time-reversal 
invariant systems and, furthermore, to those time-reversal in- 
variant systems which satisfy T^ — —1, where T is the anti- 
unitary operator generating time-reversal. Then, following 
Ref. 27, we write TuiT^^ = {Ji)ijaj. The matrix Ji is 



antisymmetric and satisfies J^ 
Hamiltonian requires 



1. T-invariance of the 



JiA = -AJi 



(17) 



^ which may, consequently, be viewed 
Multplication of w G M^^ by a quaternion can be 



Since Ji is antisymmetric and satisfies jf = —1, its eigen- 
values are ±«. Therefore, Ji defines a complex structure on 
M."^^ which may, consequently, be viewed as C^. Now con- 
sider A, which is also antisymmetric and satisfies A? — —\, 
in addition to anticommuting with Ji. It defines a quater- 
nionic structure on 
as H^/2 

defined as: (a + 6z + cj + d'k)v = av + bJiv + cAv + dJiAv. 
The possible choices of A can be obtained from a fixed one 
by performing rotations of C^, modulo those rotations which 
respect the quaternionic structure, namely Sp(A^/2). Thus, the 
set of time-reversal-invariant zero-dimensional free fermionic 
Hamiltonians with T^ ~ —1 is topologically-equivalent to 
U(iV)/Sp(7V/2). Since 7ro(U(Ar)/Sp(iV/2)) = 0, any such 
Hamiltonian can be continuously deformed into any other. 
This can be understood as a result of Kramers doubling: there 
must be an even number of fermions in the ground state so the 
division into two classes of the previous case does not exist 
here. 



C. 2D Systems: T-breaklng superconductors 

Now, let us consider systems in more than zero dimensions. 
Once again, we will assume that charge is not conserved, and 
we will also assume that time-reversal symmetry is not pre- 
served. For the sake of concreteness, let us consider a single 
band of spin-polarized electrons on a two-dimensional lattice. 
Let us suppose that the electrons condense into a (fully spin- 
polarized) Pa- -wave superconductor. For fixed superconduct- 
ing order parameter, the low-energy theory is a free fermion 
Hamiltonian for gapless fermionic excitations at the nodal 
points ±kp = (0, ±pi?). We now ask the question, what 
other order parameters could develop which would fully gap 
the fermions? For fixed values of these order parameters, we 
have a free fermion Hamiltonian. Thus, these different possi- 
ble order parameters correspond to different possible gapped 
free fermion Hamiltonians. 

The low-energy Hamiltonian of a fully spin-polarized p^- 
wave superconductor can be written in the form: 



H = tp' {iv^dxTx + ivpOyTz) tp 



(18) 



of the gap near the node, 
particle-hole space: 



The Pauli matrices t act in the 



\ -kp+k 



(19) 



This Hamiltonian is invariant under the U(l): if) -^ e ip 
which corresponds to conservation of momentum in the py 
direction (not to charge conservation). Since we will be con- 
sidering perturbations which do not respect this symmetry, it 
is convenient to introduce Majorana fermions Xi7 X2 accord- 
ing to f/; ==Xi +iX2- Then 



H = iXa (vaOxTx + VpdyT^) Xa 



(20) 



with a ~ 1,2. Note that we have suppressed the particle-hole 
index on which the Pauli matrices r act. Since Xi, X2 ^^ ^^'^h 
a 2-component real spinor, this model has 4 real Majorana 
fields. 

We now consider the possible mass terms which we could 
add to make this Hamiltonian fully gapped: 

H = iXa {vaOxTx + VpdyT^) Xa + ^Xa^abXb (21) 

If we consider the possible order parameters which could de- 
velop in this system, it is clear that there are only two choices: 
an imaginary superconducting order parameter ipy (which 
breaks time-reversal symmetry) and charge density-wave or- 
der (CDW). These take the form: 



Kr = A 



iPy 



r^Sab 



and 



M, 



CDW 
ab 



= P2kp ^y (cos e fii^ + sin e ^i^,) 



(22) 



(23) 



where vp, f a are, respectively, the Fermi velocity and slope 



where /i^'^ are Pauli matrices and 9 is an arbitrary angle. For 
an analysis of the possible mass terms in the more complex 
situation of graphene-like systems, see, for instance, Ref 32. 

Let us consider the space of mass terms with a fixed energy 
gap A which is the same for all 4 of the Majorana fermions 
in the model (i.e. a flattened mass spectrum). An arbitrary 
gapped Hamiltonian can be continuously deformed to one 
which satisfies this condition. Then we can have Aii, = ±A, 
P2k = or P21. = A, Aip = (in the latter case, arbi- 
trary 6 is allowed). If both order parameters are present, then 
the energy gap is not the same for all fermions. It's not that 
there's anything wrong with such a Hamiltonian - indeed, one 
can imagine a system developing both kinds of order. Rather, 
it is that such a Hamiltonian can be continuously deformed to 
one with either A^p = or p^j, = without closing the 
gap. For instance, if A^p > P2i^ , then the Hamiltonian can 
be continuously deformed to one with P2f. — 0. (However 
if we try to deform it to a Hamiltonian with A^p = 0, the 
gap will close at A^p — P2f. .) Hence, we conclude that 
the space of possible mass terms is topologically-equivalent 
to the disjoint union U(1)UZ2: a single one-parameter family 
and two disjoint points. 

Since 7ro(U(l) U Z2) = Z3, there are three distinct classes 



of quadratic Hamiltonians for 4 flavors of Majorana fermions 
in 2D. The one-parameter family of CDW-ordered Hamilto- 
nians counts as a single class since they can be continuously 
deformed into each other The parameter 6 is the phase of 
the CDW, which determines whether the density is maximum 
at the sites, the midpoints of the bonds, or somewhere in be- 
tween. It is important to keep in mind, however, that, although 
there is no topological obstruction to continuously deforming 
one 6 into another, there may be an energetic penalty which 
makes it costly. For instance, the coupling of the system to 
the lattice may prefer some particular value of 6. The classi- 
fication discussed here accounts only for topological obstruc- 
tions; the possibility of energetic barriers must be analyzed by 
different methods. 

We can restate the preceding analysis in the following, more 
abstract language. This formulation will make it clear that we 
haven't overlooked a possible mass term and will generalize 
to more complicated free fermion models. Let us write ^^ = 
TxSab, 72 = Tjab- Then 



{7„7j = 2,5, 



(24) 



The Dirac Hamiltonian for iV := 4 Majorana fermion fields 
takes the form 



H = ixilrdr + M)x 



(25) 



The matrix M plays the role that the matrix A did in the zero- 
dimensional case. As in that case, we assume a flattened spec- 
trum which here means that each Majorana fermion field has 
the same gap and that this gap is equal to 1. (It does not mean 
that the energy is independent of the momentum k.) In order 
to satisfy this, we must require that 



{7^,M} = and M^ 



-1 



(26) 



Note that it is customary to write the Dirac Hamiltonian in 
a slightly different form, 

H = ip{i-fidi + ■m)'ip (27) 

which can be massaged into the form of (25) using ip = tp^jo'. 

H = V^(no7i5i + "^7o)'0 
= ^p' {iaidi + to/3)')/' 
= itp^aidi — im(3)ip (28) 



where ai ~ 797^ and /3 = 70. Thus, if we write 7^ = ai and 
M = —imj3 and consider Majorana fermions (or decompose 
Dirac fermions into Majoranas), we recover (25). We have 
used the form (25) so that it is analogous to (13), with (7^9^ + 
M) replacing Aij and the i pulled out front. Then, the matrix 
M determines the gaps of the various modes in the same way 
as A does in the zero-dimensional case. Similarly, assuming a 
'flattened' spectrum leads to the condition M^ = — 1. 

How many ways can we choose such an Ml Since 7I = 1, 
its eigenvalues are ±1. Hence, viewed as a linear map from 
M** to itself, this matrix divides W^ into two 2D subspaces 
R** = X_|_ © X^ with eigenvalue ±1 under 73, respectively. 



For 72 = TzSab, this is trivial: 

l\ ^. f 1 \ (I 
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/ ' V 



(29) 



where Tz acts on the first spinor and the second spinor is in- 
dexed by a = 1, 2, i.e. is acted on by the Pauli matrices /i^'^ 
in (23). This construction generalizes straightforwardly to ar- 
bitrary numbers A^ of Majorana fermions, which is why we 
use it now. 

Now 7i-M commutes with 73 and satisfies (7]^M)^ = 1. 
Thus, it maps X+ to itself and defines subspaces X\_ , X^ with 
eigenvalue ±1 under ^^Ad (and equivalently for XJ). X^ 
can decomposed into X\_ ffi X"^ = X+. Choosing M is thus 
equivalent to choosing a linear subspace X\_ of X+. 

This can be divided into three cases. If JiM has one posi- 
tive eigenvalue and one negative one when acting on X+ then 
the space of possible choices of JiM is equal to the space of 
ID linear subspaces of M^, which is simply U(l). If, on the 
other hand, 7]^M has two positive eigenvalues, then there is 
a unique choice, which is simply M = 75^73. If 7]^M has 
two negative eigenvalues, then there is again a unique choice, 
M — — 7i72- Therefore, the space of possible Ms is topolog- 
ically equivalent to U(l) U Z2. 

Now, suppose that we have 2A^ Majorana fermions. Then 
72 defines iV-dimensional eigenspaces X+ , X^ such that 
X_(_ © X_ and JiM defines eigenspaces of X^: 
X+. If 7iM has k positive eigenvalues and 
N ~ k negative ones, then the space of possible choices of 
7jM is 0(N)/0(k)xO(N-k), i.e we can take the restriction of 
7j^ M to X^ to be of the form 



p2JV 



x\ 



x\ 



/I 



7^/Vf = O^ 
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O (30) 



-1/ 



with k diagonal entries equal to +1 and N — k entries equal 
to —1. Thus, the space of Hamiltonians for N flavors of free 
Majorana fermions is topologically equivalent to 



N 



M 



2N 



\JO{N)/{0{k) xO{N-k)) 



(31) 



k=0 



However, since 7ro(0(/V)/(0(fc) x 0(/V- fc))) = 0, indepen- 
dent of k (note that is the group with a single element, not 
the empty set 0), TToiM2N) — ^n+i- 

In the model analyzed above, we had only a single spin- 
polarized band of electrons. By increasing the number of 
bands and allowing both spins, we can increase the number 
of flavors of Majorana fermions. In principle, the number of 
bands in a solid is infinity. Usually, we can introduce a cutoff 
and restrict attention to a few bands near the Fermi energy. 
However, for a purely topological classification, we can ig- 
nore energetics and consider all bands on equal footing. Then 
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we can take N -^ oo, so that 7ro(A^co) = Z. This classifi- 
cation permits us to deform Hamiltonians into each other so 
long as there is no topological obstruction, with no regard to 
how energetically costly the deformation may be. Thus, the 
2N — 4 classification which we discussed above can perhaps 
be viewed as a 'hybrid' classification which looks for topolog- 
ical obstructions in a class of models with a fixed set of bands 
close to the Fermi energy. 

But even this point of view is not really natural. The dis- 
cussion above took as its starting point an expansion about a 
Px superconductor; the p^. superconducting order parameter 
was assumed to be large and fixed while the ipy and CDW 
order parameters were assumed to be small. In other words, 
we assumed that the system was at a point in parameter space 
at which the gap, though non-zero, was small at two points in 
the Brillouin zone (the intersection points of the nodal line in 
the px superconducting order parameter with the Fermi sur- 
face). This allowed us to expand the Hamiltonian about these 
points in the Brillouin zone and write it in Dirac form. And 
this may, indeed, be reasonable in a system in which px super- 
conducting order is strong (i.e. highly energetically-favored) 
and other orders are weak. However, a topological classifica- 
tion should allow us to take the system into regimes in which 
Px superconductivity is small and other orders are large. Sup- 
pose, for instance, that we took our model of spin-polarized 
electrons (which we assume, for simplicity, to be at half -filling 
on the square lattice) and went into a regime in which there 
was a large (i2.2_j,2 -density-wave (or 'staggered flux') order 

parameter-'-' (cj^ioCk) = i^{coskxa — coskyO), where a 
is the lattice constant and $ is the magnitude of the order 
parameter. With nearest-neighbor hopping only, the energy 
spectrum is £'^ = (2t)^(cosfca;0 + cosfcya)^ + <I>^(cosfca;a — 
cos fej,o)^. Thus, the gap vanishes at 4 points, (±7r/2, ±7r/2) 
and (=F7r/2, ±7r/2). The Hamiltonian can be linearized in the 
vicinity of these points: 

H = tpl {ivAdxTx + ivpdyTz) V'l 

+ 'ipl {iv\dyTx + ivpdxTz) ip2 (32) 

where vp, va are, respectively, the Fermi velocity and slope 
of the gap near the nodes; the subscripts 1,2 refer to the two 
sets of nodes (±7r/2, ±7r/2) and (=F7r/2,±7r/2); and tpA, 
A = 1,2 are defined by: 



Adk) 



^iTr/2,±Tr/2)+k 
"i-7T/2,TTT/2)+k 



(33) 



If we introduce Majorana fermions ipA — Xai + *Xa2' ^^^^ 
we can write this Hamiltonian with possible mass terms as: 

H = iXla {vaOxTx + VpdyT^) Xla 

+ iX2a {vAdyTx + VpdxT^) X2a 

+ ^XAa MAa,Bb Xsb (34) 

We have suppressed the spinor indices (e.g. xn is a two- 
component spinor); with these indices included, MAa,Bb is 
an 8 X 8 matrix. However, in order for the gap to be the 



same for all flavors, the mass matrix must anticommute with 
Tx,z- Thus, MAa,Bb = TyMAaMb- The matrix MAa,Bb 
can have 0, 1,2,3, or 4 eigenvalues equal to +1 (with the 
rest being —1). The spaces of such mass terms are, re- 
spectively, 0, 0(4)/(0(l) X 0(3)), 0(4)/(0(2) x 0(2)), 
0(4)/(0(3) X 0(1)), and 0. Mass terms with or 4 eigen- 
values equal to +1 correspond physically to ztidxy-density 

wave order, (cj^+qCk) = ±sinfca;a sinfcj,a. Mass terms with 
2 eigenvalues equal to +1 correspond physically to supercon- 
ductivity, to Q' — (tt, 0) CDW order, and to linear combina- 
tions of the two. Mass terms with 1 or 3 eigenvalues equal 
to +1 correspond to (physically unlikely) hybrid orders with, 
for instance, superconductivity at (±7r/2, ±7r/2) and ztidxy- 
density wave order at (±7r/2, =F7r/2). Clearly, this is the 
2N = 8 case of the general classification discussed above. 
Thus, the same underlying physical degrees of freedom - a 
single band of spin-polarized electrons on a square lattice - 
can correspond to either 2N = 4 or 2N — 8, depending on 
where the system is in parameter space. One can imagine re- 
gions of parameter space where the gap is small at an arbitrary 
number A^ of points. Thus, if we restrict ourselves to systems 
with a single band, then different regions of the parameter 
space (with different numbers of points at which the gap is 
small) will have very different topologies. Although such a 
classification may be a necessary evil in some contexts, it is 
far preferable, given the choice, to allow topology to work un- 
fettered by energetics. Then, we can consider a large number 
n of bands. Suppose that the gap becomes small at r points in 
the Brillouin zone in each band. Then, the low-energy Hamil- 
tonian takes the Dirac form for 2A^ = 2rn Majorana fermion 
fields. As we will see below, if N is sufficiently large, the 
topology of the space of possible mass terms will be indepen- 
dent of N. Consequently, for n sufficiently large, the topology 
of the space of possible mass terms will be independent of r. 
In other words, we are in the happy situation in which the 
topology of the space of Hamiltonians will be the same in the 
vicinity of any gap closing. But any gapped Hamiltonian can 
be continuously deformed so that the gap becomes small at 
some points in the Brillouin zone. Thus, the problem of clas- 
sifying gapped free fermion Hamiltonians in rf-dimensions is 
equivalent to the problem of classifying possible mass terms 
in a generic d-dimensional Dirac Hamiltonian so long as the 
number of bands is sufficiently large^^. This statement can be 
made more precise and put on more solid mathematical foot- 
ing using ideas which we discuss in Appendix B. 



D. Classification of Topological Defects 

The topological classification described above holds not 
only for classes of translationally-invariant Hamiltonians such 
as (25), but also for topological defects within a class. Sup- 
pose, for instance, that we consider (25) with a mass which 
varies slowly as the origin is encircled at a great distance. We 
can ask whether such a Hamiltonian can be continuously de- 
formed into a uniform one. In a system in which the mass term 
is understood as arising as a result of some kind of underly- 
ing ordering such as superconductivity or CDW order, we are 
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simply talking about topological defects in an ordered media, 
but with the caveat that the order parameter is allowed to ex- 
plore a very large space which may include many physically 
distinct or unnatural orders, subject only to the condition that 
the gap not close. 

Let us, for the sake of concreteness, assume that we have 
a mass term with N/2 positive eigenvalues when restricted 
to the +1 eigenspace of 72. (For N large, the answer obvi- 
ously cannot depend on the number of positive eigenvalues k 
so long as k scales with N . Thus, we will denote the space 
M2N defined in Eq. 31 by Z x 0{N)/{0{N/2) x 0{N/2)) 
where the integers in Z correspond to the number of pos- 
itive eigenvalues of the mass term when restricted to the 
+1 eigenspace of 72.) Then M{r — 00, 9) defines a loop 
in 0{N)/{0{N/2) x 0{N/2)) which cannot be continu- 
ously unwound if it corresponds to a nontrivial element of 
^i(0(iV)/(0(iV/2)xO(A^/2))). 

To compute 7ri(0(iV)/(0(iV/2) x 0(A^/2))), we 
parametrize 0{N)/{0{N/2) x 0{N/2)) by symmetric 
matrices K which satisfy K^ = 1 and tr(if) = 0. (Such 
matrices decompose R^ into their +1 and —1 eigenspaces: 
R^ :=:V+®V-. iC can be written in the form: K = O^KaO, 
where Kq has N/2 diagonal entries equal to +1 and N /2 
equal to -1, i.e Kq = diag(l, . . . , 1, -1, . . . , -1)-) Note that 
any such K is itself an orthogonal matrix, i.e. an element of 
0(iV); thus 0(iV)/(0(iV/2) x 0(7V/2)) can be viewed as a 
submanifold of 0(iV) in a canonical way. Consider a curve 
L(A) in 0(A^)/(0(A^/2) x 0(iV/2)) with L(0) = K and 
L(7r) — —K. We will parametrize it as L{\) = K e^^, where 
A is in the Lie algebra of O(A^). In order for this loop to re- 
main in 0{N)/{0{N/2) X 0(iV/2)), we need {K e^^f = 1. 
Since [K e^'^f = Ke^^Ke'^^ = ^iXKAK ^XA^ ^j^j^ 
condition implies that KA = —AK. In order to have 
L(7r) ~ —K, we need A^ = —1. Such a curve is, in fact, a 
minimal geodesic from K to —K. Each such geodesic can 
be represented by its midpoint L{'k/2) = KA, so the space 
of such geodesies is equivalent to the space of matrices A 
satisfying A^ — —1 and KA = ~AK. As discussed in 
Ref 28, the space of minimal geodesies is a good enough 
approximation to the entire space of loops (essentially 
because an arbitrary loop can be deformed to a geodesic) 
that we can compute tti from the space of minimal geodesies 
just as well as from the space of loops. Thus, the loop 
space of 0{N/2)/{0{k) x 0{N/2 - k)) is homotopically 
equivalent to the space of matrices A satisfying A^ = —1 
and KA = —AK. Since it anticommutes with K, KA maps 
the +1 eigenspace of K to the —1 eigenspace. It is clearly a 
length-preserving map since (KA)'^ ~ 1 and, since the ±1 
eigenspaces of K are isomorphic to M^/^, KA defines an ele- 
ment of 0(A^/2). ThusaloopinO(A^)/(0(A^/2) x 0(iV/2)) 
corresponds to an element of 0(A^/2) or, in other words: 

7r,{0{N)/{0{N/2) X 0{N/2))) - 7ro(0(iV/2)). (35) 

The latter group is simply Z2 since 0{N/2) has two con- 
nected components: (1) pure rotations and (2) rotations com- 
bined with a reflection. 

It might come as a surprise that we find a Z2 classification 



for point-like defects in two dimensions. Indeed, if we require 
that the superconducting order parameter has fixed amplitude 
at infinity, then vortices of arbitrary winding number are stable 
and we have a Z classification. However, in the classification 
discussed here, we require a weaker condition be satisfied: 
that the fermionic gap remain constant. Consequently, a vor- 
tex configuration of even winding number can be unwound 
without closing the free fermion gap by, for instance, 'rotat- 
ing' superconductivity into CDW order. 



E. 3D Systems with No Symmetry 



With these examples under our belts, we now turn to the 
case which is of primary interest in this paper: free fermion 
systems in three dimensions without time-reversal or charge- 
conservation symmetry. We consider the Dirac Hamiltonian 
in 31? for an 2A^-component Majorana fermion field %: 



H = ix. (9i7i + ^272 + ^373) X + ixMx 



(36) 



In the previous section, we discussed a lattice model which 
realizes (36) in its continuum limit with 2N = 8. Different 
mass terms correspond to different quadratic perturbations of 
this model which open a gap (which can be viewed as order 
parameters which we are turning on at the mean-field level). 
We could classify such terms by considering, from a physical 
perspective, all such ways of opening a gap. However, we 
will instead determine the topology of the space of mass terms 
(and, thereby, the space of gapped free fermion Hamiltonians) 
by the same mathematical methods by which we analyzed the 
2D case. 



Since 7^ = 1 and has vanishing trace, this matrix decom- 
poses M^^ into its ±1 eigenspaces: M^^ = X+ ® X_. Now 
(7273)^ = -1 and [71,7273] == 0- Therefore, 7273 is a com- 
plex structure on X+ (and also on X_), i.e. we can define 
multiplication of vectors v <E X^ by complex scalars accord- 
ing to (a + bi)v = av + ^2l?,'^- (Consequently, we can view 
X+ as C^/^.) Now, consider a possible mass term M, with 
IVP = -1. (72M)2 = 1 and [71,72^] = 0. Let Y be 
the subspace of X^ with eigenvalue +1 under 72M. Since 
{72M, 7273} = 0, 7273^ is the subspace of X+ with eigen- 
value — 1 under 72M . In other words, X+ — Y ® l-fliX ^ 
i.e. y is a real subspace of X+. Hence, the space of choices 
of M is the space of real subspaces Y C X^ (or, equiva- 
lently, of real subspaces M^/^ c C^/^). Given any fixed real 
subspace Y C X+, we can obtain all others by performing 
U(Af/2) rotations of X+, but two such rotations give the same 
real subspace if they differ only by an 0{N/2) rotation of Y . 
Thus, the space of gapped Hamiltonians for 2N free Majo- 
rana fermion fields in iD with no symmetry is topologically- 
equivalent to U(iV/2)/0(A^/2). In the remaining sections of 
this paper, we will be discussing topological defects in such 
systems and their motions. 
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F. General Classification and Bott Periodicity 

Before doing so, we pause for a minute to consider the clas- 
sification in other dimensions and in the presence of symme- 
tries such as time-reversal and charge conservation. We have 
seen that systems with no symmetry in d = 0, 2, 3 are clas- 
sified by the spaces 0(2iV)/U(A^), Z x o(jv/°/:^V/2) ' ^'^ 
U(iV/2)/0(7V/2). By similar methods, it can be shown that 
the d — 1 case is classified by O(A^). As we have seen, in- 
creasing the spatial dimension increases the number of 7 ma- 
trices by one. The problem of choosing 71 , . . . , 7^ satisfying 
{lijlj} = 2(5ij and M which anti-commutes with the 7iS and 
squares to —1 leads us to subspaces of M^^ of smaller and 
smaller dimension, isometrics between these spaces, or com- 
plex of quatemionic structures on these spaces. This leads the 
progression of spaces in the top row of Table 1. 

At the same time, we have seen that a time-reversal- 
invariant system in d = is classified by U(7V)/Sp(7V/2). 
Suppose that we add a discrete anti-unitary symmetry 
SiS^^ = —i defined by 



SaiS ^ = {J)ijaj 



(37) 



which squares to J^ = —1. It must anti-commute with the 
mass term 



{J, M} = 



(38) 



in order to ensure invariance under the symmetry, so choos- 
ing a J amounts to adding a complex structure, which leads 
to the opposite progression of classifying spaces. Consider, 
as an example of the preceding statements, a time-reversal in- 
variant system in d = 3. Then time-reversal symmetry T is an 
example of a symmetry generator J discussed in the previous 
paragraph. We define a real subspace Y C X+, in a similar 
manner as above, but now as the subspace of X+ with eigen- 
value +1 under 72T, rather than under 72M. Once again, 
X+ = r 07273^- Now, {73^,72^} = 0, and {^.M)^ = 1, 
so the +1 eigenspace of 73M is a linear subspace of Y. The 



set of all such linear subspaces is Z x 



0(JV/2) 
0(Ar/4)xO(Ar/4) 



But this 



is the same classifying space as for a system with no symme- 
try in rf = 2 (apart from a reduction of iV by a factor of 2). 
Thus, we are led to the list of classifying spaces for gapped 
free fermion Hamiltonians in Table 1. 

In this table, Q refers to charge-conservation symme- 
try. Charge conservation is due to the invariance of the 
Hamiltonian of a system under the U(l) symmetry Ci -^ 
e^^Ci. In terms of Majorana fermions a^ defined according 
to Cj = (a2j-i + *a2j)/2, the symmetry takes the form 
a2j-i — > cos9a2j-i + sm6a2j, a2j — > — siii9a2j-i + 
cos 0a2j ■ However, if a free fermion Hamiltonian is invari- 
ant under the discrete symmetry q — > ici or, equivalently, 
a2j-i ^ 0'2j, 0-2] ^ ^o,2j-i7 then it is automatically in- 
variant under the full U(l) as well, and conserves charge^^. 
Thus, we can treat charge conservation as a discrete symme- 
try Q which is unitary, squares to —1, and commutes with the 
Hamiltonian (i.e. with the 7 matrices and M). Since Q trans- 
forms Ci — > iCi, it anti-commutes with T. Note further that 



if a system has time-reversal symmetry, then the product of 
time-reversal T and charge conservation Q is a discrete anti- 
unitary symmetry, QT which anti-commutes with the Hamil- 
tonian and with T and squares to —1. Then QT is defined by 
a choice of matrix J, analogous to T, as in Eq. 37. If the sys- 
tem is not time-reversal-invariant, then charge conservation is 
a unitary symmetry. It is easier then to work with complex 
fermions, and the classification of such systems falls into an 
entirely different sequence, as discussed in Appendix A.) 

If a system is both time-reversal symmetric and charge- 
conserving, i.e. if it is a time-reversal invariant insulator, then 
it may have an additional symmetry which guarantees that the 
eigenvalues of the Hamiltonian come in ±E pairs, just as in a 
superconductor. An example of such a symmetry is the sublat- 
tice symmetry of Hamiltonians on a bipartite lattice in which 
fermions can hop directly from the A sublattice to the B sub- 
lattice but cannot hop directly between sites on the same sub- 
lattice. In such a case, the system is invariant under a unitary 
symmetry x defined as follows. If we block diagonalize x so 
that one block acts on sites in the A sublattice and the other on 
sites in the B sublattice, then we can write x = diag(fc, — fc), 
i.e. ai(x) — > —kijaj{x.) for x e A and ai{x.) — > kijaj{x.) for 
X G B. This symmetry transforms the Hamiltonian to minus 
itself if k"^ — \ or, in other words, if x^ = 1. Then xQ is a uni- 
tary symmetry which squares to —1 and anti-commutes with 
the Hamiltonian, T, and QT. Hence xQ, too, is defined by a 
choice of matrix J, as in Eq. 37. We will call such a symmetry 
a sublattice symmetry x and a system satisfying this symme- 
try a 'bipartite' system, but the symmetry may have a different 
microscopic origin. 

In an electron system, time-reversal ordinarily squares to 
— 1, because the transformation law is c-|- -^ c^, c^ ^- ~c^, 
as we have thus far assumed in taking J^ — —1. However, it 
is possible to have a system of fully spin-polarized electrons 
which has an anti -unitary symmetry T which squares to +1. 
(One might object to calling this symmetry time-reversal be- 
cause it doesn't reverse the electron spins, but T is a natural 
label because it is a symmetry which is just as good for the 
present purposes.) Then, since J^ = 1, a choice of J is sim- 
ilar to a choice of a 7 matrix. In general, symmetries (37) 
which square to +1 have the same effect on the topology of 
the space of free fermion Hamiltonians as adding dimensions 
since each such J defines a subspace of half the dimension 
within the eigenspaces of the 7 matrices. This is true for sys- 
tems with T^ = 1. 

SU(2) spin-rotation-invariant and time-reversal-invariant 
insulators (systems with T and Q) effectively fall in this cate- 
gory. The Hamiltonian for such a system can be written in the 
form H = h(E)l2 where the second factor is the 2x2 identity 
matrix acting on the spin index. Then time-reversal can be 
written in the form T = t (E) ia-y, where i^ = 1, and Q can be 
written in the form Q — q(E)l2,S'0 that QT = qt(E) ia-y, where 
(gt)^ = 1. Thus, since the matrix i<Ty squares to —1, the sym- 
metries T and QT have effectively become symmetries which 
square to +1. They now move the system through the progres- 
sion of classifying spaces in the same direction as increasing 
the dimension, i.e. in the opposite direction to symmetries 
which square to —1. Thus, SU(2) spin-rotation-invariant and 
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dim.: 

SU(2), T, Q 

SU(2), T, Q, X 

no symm. 

Tonly 

TandQ 

T,Q,x 



Z X 



J) 

0(N) 



0{N/2)xO{N/2) 
0(iV/4) 

0(2A^)/U(A^) 
U(7V)/Sp(iV/2) 

Sp(JV/2) 
Sp(Af/4)xSp(Af/4) 

Sp(7V/4) 



Z X 



1 

U(iV/2)/0(iV/2) 

^ 0(JV/4) 

^ ^ 0(Af/8)xO(Af/8) 

0(JV) 

0(iV)/U(iV/2) 

U(iV/2)/Sp(iV/4) 



Sp(iV/4)/U(iV/4) 

U(iV/8)/0(iV/8) 

Z X 



0(JV) 



Z X 



Sp(iV/4) 
Sp(Ar/8)xSp(Ar/8) 



0(Ar/2)xO(Ar/2) 
0{N/2) 

0(iV/2)/U(iV/4) 
U(iV/4)/Sp(iV/8) 



Sp(7V/8) 

Sp(7V/16)/U(iV/16) 

U(iV/2)/0(7V/2) 

7 X 0(^/2) 

^ ^ 0(Af/4)xO(7V/4) 

0(iV/4) 

0(Ar/4)/U(7V/8) 



Z X 



Sp{N/S) 



Sp(JV/16)xSp(Ar/16) 

Sp(iV/32) 
Sp(iV/4)/U(iV/4) 



TABLE I: The period-8 (in both dimension and number of symmetries) table of classifying spaces for free fermion Hamiltonians for N complex 
— 2N real (Majorana) fermion fields in dimensions d = 0, 1, 2, 3, . . . with no symmetries; time-reversal symmetry (T) only; time-reversal 
and charge conservation symmetries (T and Q); time-reversal, charge conservation, and sublattice symmetries (T, Q, and x)\ and the latter 
two cases with SU(2) symmetry. As a result of the period-8 nature of the table, the top two rows could equally well be the bottom two rows of 
the table. Moving p steps to the right and p steps down leads to the same classifying space (but for 1/2'' as many fermion fields), which is a 
reflection of Bott periodicity, as explained in the text. The number of disconnected components of any such classifying space - i.e. the number 
of different phases in that symmetry class and dimension - is given by the corresponding tto, which may be found in Eq. 40. Higher homotopy 
groups, which classify defects, can be computed using Eq. 39. Table III, given in Appendix A, is the analogous table for charge-conserving 
Hamiltonians without time-reversal symmetry. 



time-reversal-invariant insulators in d dimensions are classi- 
fied by the same space as systems with no symmetry in d + 2 
dimensions. However, in a system which, in addition, has sub- 
lattice symmetry x = x (E) I2, "^^ have (qx)^^. Thus, sublat- 
tice symmetry is still a symmetry which squares to —1. Since 
the two symmetries which square to +1 (T and QT) have the 
same effect as increasing the dimension while the symmetry 
which squares to —1 has the same effect as decreasing the 
dimension, SU(2) spin-rotation-invariant and time-reversal- 
invariant insulators with sublattice symmetry in d dimensions 
are classified by the same space as systems with no symmetry 
in d + 1 dimensions (but with N replaced by N/4). Simi- 
lar considerations apply to superconductors with SU(2) spin- 
rotational symmetry. 

In order to discuss topological defects in the sys- 
tems discussed here, it is useful to return to the ar- 
guments which led to (35). By showing that the 
space of loops in 0(7V/2)/(0(A^/4) x 0{N/4)) is 
well-approximated by 0(A^/4), we not only showed 
that 7ri(0(A^/2)/(0(7V/4) x 0{N/4:))) = 7ro(0(iV/4)) 
but, in fact, that 7rfe(0(A^/2)/(0(7V/4) x 0(iV/4))) == 
7rfc_i(0(iV/4)) (see Ref. 28). Continuing in the same 
way, we can approximate the loop space of 0(A^/4) (i.e. 
the space of loops in 0(iV/4)) by minimal geodesies from 
I to -I: L'{X) = e^-^i where Aj = -1. The mid- 
point of such a geodesic, L'{tt/2) — A\ again defines 
a complex structure A^ ~ O^ JO, where J is given by 
(16) so that 7rfe(0(iV/4)) == TTk-i{0{N/4)/\J{N/%)). In 
a similar way, minimal geodesies in 0{N /A)/\]{N /^) from 
Ai to —Ai can be parametrized by their mid-points A2, 
which square to —1 and anti-commute with Ai, thereby 
defining a quatemionic structure, so that the loop space of 
0(7V/4)/U(7V/8) is equivalent to U(iV/8)/Sp(A^/8) and, 
hence 7rfc(0(7V/4)/U(7V/8)) = ^fe_i(U(iV/8)/Sp(iV/8)). 
Thus, we see that the passage from one of the classifying 
spaces to its loop space is the same as the imposition of a 



symmetry such as time-reversal to a system classified by that 
space: both involve the choice of successive anticommuting 
complex structures. Continuing in this fashion (see Ref. 28), 
we recover Bott periodicity: 

7rfe(0(167V)) ^ 

^fe_i(0(167V)/U(8iV)) = 7rfc_2(U(87V)/Sp(4iV)) 

= ^fe_3(ZxSp(4iV)/(Sp(27V)xSp(2Af))) = ^fc_4(Sp(2A^)) 

= ^fe_5(Sp(2iV)/U(2Ar)) = ^fc_6(U(2Af)/0(2iV)) 

= 7rfe_7(Z X 0{2N)I{0{N) x 0{N))) 

- 7rfe_8(0(7V)) (39) 

The approximations made at each step require that N be in 
the stable limit, in which the desired homotopy groups are 
independent of N . For instance, TTk{0{N)) is independent of 

iV for TV > k/2. 

It is straightforward to compute ttq for each of these groups: 



^o(0(iV)) = Z2 

7ro(0(2iV)/U(7V)) = Z2 

7ro(U(2iV)/Sp(7V)) = 

7ro(Z X Sp(27V)/Sp(iV) x Sp(iV)) = Z 

7ro(Sp(iV)) = 

7ro(Sp(iV)/U(iV)) = 

7ro(U(iV)/0(iV)) - 

7ro(Z X 0(2iV)/(0(iV) x 0(iV))) = Z 



(40) 



Combining (40) with (39), we can compute any of the stable 
homotopy groups of the above 8 classifying spaces. As dis- 
cussed above, the space of gapped free fermion Hamiltonians 
in d-dimensions in a given symmetry class (determined by the 
number modulo 8 of symmetries squaring to —1 minus the 
number of those squaring to +1) is homotopically-equivalent 
to one of these classifying spaces. Thus, using (40) with 
(39) to compute the stable homotopy groups of these classi- 
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Symmetry classes 


Physical realizations 


d = l 


d = 2 


d = 3 


D 


SC 


p-wave SC 


{p + ip)-SC 





Dill 


TRISC 


Z2 


{p + ip){p- ip)-SC 


He3-B 


All 


TRl ins. 





HgTe Quantum well 


Bii_xSbx, Bi2Sc3, etc. 


Cll 


Bipartite TRl ins. 


Carbon nanotube 





Z2 


C 


Singlet SC 





{d + id)-SC 





CI 


Singlet TRl SC 








z 


Al 


TRl ins. w/o SOC 











BDl 


Bipartite TRl ins. w/o SOC 


Carbon nanotube 









TABLE II: Topological periodic table in physical dimensions 1, 2, 3. The first column contains 8 of the 10 symmetry classes in the Cartan 
notation adopted by Schnyder et al?^ , following Zirnbauer'*"^^ . The second column contains the requirements for physical systems which 
can realize the corresponding symmetry classes, "w/o" stands for "without". SC stands for superconductivity, TRl for time-reversal invariant, 
and SOC for spin-orbit coupling. The three columns d = 1, 2, 3 list topological states in the spatial dimensions 1, 2, 3 respectively. means 
the topological classification is trivial. The red labels Z and Z2 stand for topological states classified by these groups but for which states 
corresponding to non- trivial elements of Z or Z2 have not been realized in realistic materials. The blue text stands for topological states for 
which a well-defined physical model has been proposed but convincing experimental candidate has not been found yet. (See text for more 
discussions on the realistic materials.) 



fying spaces leads to a complete classification of topological 
states and topological defects in all dimensions and symmetry 
classes, as we now discuss. 

Gapped Hamiltonians with a given symmetry and dimen- 
sion are classified by ttq of the corresponding classifying 
space in Table 1. Due to Bott periodicity, the table is pe- 
riodic along both directions of dimension and symmetry, so 
that there are 8 distinct symmetry classes. Ryu et. af^ de- 
noted these classes using the Cartan classification of symmet- 
ric spaces, following the corresponding classification of dis- 
ordered systems and random matrix theory^^'^^ which was ap- 
plied to the (potentially-gapless) surface states of these sys- 
tems. In this notation, systems with no symmetry are in class 
D, those with T only are in Dill, and those with T and Q 
are in All. The other 5 symmetry classes, C, CI, Cll, Al, and 
BDl arise, arise in systems which have spin-rotational sym- 
metry or a sublattice symmetry. There are actually 2 more 
symmetry classes (denoted by A and Alll in the random ma- 
trix theory) which lie on a separate 2x2 periodic table, which 
is less relevant to the present work and will be discussed in 
the Appendix A. In Table 11 we have listed examples of 
topologically-nontrivial states in physical dimensions 1,2,3 in 
all 8 symmetry classes. To help with the physical understand- 
ing of these symmetry classes, we have also listed the phys- 
ical requirements for the realization of each symmetry class. 
In each dimension, there are two symmetry classes in which 
the topological states are classified by integer invariants and 
two symmetry classes in which the different states are distin- 
guished by Z2 invariants. In all the cases in which a real ma- 
terial or a well-defined physical model system is known with 
non-trivial Z or Z2 invariant, we have listed a typical example 
in the table. In some of the symmetry classes, non-trivial ex- 
amples have not been realized yet, in which case we leave the 
topological classification Z or Z2 in the corresponding posi- 
tion in the table. 

In one dimension, generic superconductors (class D) are 
classified by Z2, of which the nontrivial example is a p- 



wave superconductor with a single Majorana zero mode on 
the edge. The time-reversal invariant superconductors (class 
Dill) are also classified by Z2. The nontrivial example is 
a superconductor in which spin up electrons pair into a p- 
wave superconductor and spin down electron form another 
p-wave superconductor which is exactly the time-reversal of 
the spin-up one. Such a superconductor has two Majorana 
zero modes on the edge which form a Kramers pair and are 
topologically protected. The two integer classes are bipar- 
tite time-reversal invariant insulators with (Cll) and without 
(BDl) spin-orbit coupling. An example of the BDl class is 
a graphene ribbon, or equivalently a carbon nanotube with a 
zigzag edge.^^'^^. The low-energy band structure of graphene 
and carbon nanotubes is well-described by a tight-binding 
model with nearest-neighbor hopping on a honeycomb lat- 
tice, which is bipaitite. The integer-valued topological quan- 
tum number coiTesponds to the number of zero modes on the 
edge, which depends on the orientation of the nanotube. Be- 
cause carbon has negligible spin-orbit coupling, to a good ap- 
proximation it can be viewed as a system in the BDl class, 
but it can also be considered as a system in class Cll when 
spin-orbit coupling is taken into account. In two dimensions, 
generic superconductors (class D) are classified by an integer, 
coiTesponding to the number of chiral Majorana edge states 
on the edge. The first nontrivial example was the p + ip 
wave superconductor, shown by Read and Green'^ to have 
one chiral Majorana edge state. Non-trivial superconductors 
in symmetry class D are examples of topological supercon- 
ductors. Some topological superconductors can be consistent 
with spin rotation symmetry; singlet superconductors (class 
C) are also classified by integer, with the simplest physical 
example a d + id wave superconductor Similar to the ID 
case, the time-reversal invariant superconductors (class Dill) 
are classified by Z2, of which the nontrivial example is a su- 
perconductor with p + ip pairing of spin-up electrons and 
p—ip pairing of spin-down electrons. ^^■^**-^^ The other symme- 
try class in 2D with a Z2 classification is composed of time- 
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reversal invariant insulators (class All), also known as quan- 
tum spin Hall insulators""' '*^. The quantum spin Hall insula- 
tor phase has been theoretically predicted'*^^ and experimen- 
tally realized'*'* in HgTe quantum wells. In three dimensions, 
time-reversal invariant insulators (class All) are also classi- 
fied by ^2.^'"^^ The Z^ topological invariant corresponds to 
a topological magneto-electric response with quantized coef- 
ficient (? = 0, TT^"*. Several nontrivial topological insulators 
in this class have been theoretically predicted and experimen- 
tally realized, including Bii_xSbx alloy"*^'*^ and the family 
of BiaScs, BiaTca, SbaTea'*''^^ In 3D, time-reversal invari- 
ant superconductors (class Dill) are classified by an integer, 
corresponding to the number of massless Majorana cones on 
the surface. ^^ A nontrivial example with topological quantum 
number A^ = 1 turns out to be the B phase of He^.^^'-'^'^" 
The other classes with nontrivial topological classification in 
3D are singlet time-reversal invariant superconductors (CI), 
classified by an integer; and bipartite time-reversal invariant 
insulators (CII), classified by Z^- Some models have been 
proposed^' but no realistic material proposal or experimental 
realization has been found in these two classes. We would like 
to note that different physical systems can correspond to the 
same symmetry class. For example, bipartite superconductors 
are also classified by the BDI class. 

The two remaining symmetry classes (unitary (A) and chi- 
ral unitary (AIII)) corresponds to systems with charge conser- 
vation symmetry but without time-reversal symmetry, which 
forms a separate 2x2 periodic table. For the sake of com- 
pleteness, we carry out the preceding analysis for these two 
classes in Appendix A. 

Topological defects in these states are classified by higher 
homotopy groups of the classifying spaces. Following 
the convention of Ref.^^, we name the classifying spaces 
by i?g,g = 0,1, 2,..., 7, with R^ = 0(iV), i?2 = 
0(2iV)/C/(iV), ... R^ = [/(iV)/0(iV), i?o = Z x 
0(2iV)/(0(iV)xO(A^))intheorderofEq. (40). The symme- 
tries in Table II can be labeled by p = 0, 1, 2, ..., 7, so that in 
d dimensions and p-th symmetry class, the classifying space 
is R2+p-d- A topological defect with dimension D {D < d) 
is classified by 



in Ref 20, which we discussed in the strong-coupling limit in 
the previous section, is an example of this particular 'small N' 
case, which is why the hedgehogs in that model are classified 
by a winding number e Z. 



IV. EXCHANGING PARTICLES CONNECTED BY 
RIBBONS IN 3D 

In this section, we will take an 0(3) non-linear a model, 
i.e. one with target space S'^, as a toy model for our problem. 
It is essentially the '8-band' model discussed in Ref. 20 and 
in the strong-coupling limit in Section II. As noted above, it 
is the 2A^ = 8 limit of the classification reviewed above. It 
will be more familiar to most readers and easier to visualize 
than the full problem which we discuss in the next section. We 
will give a heuristic explanation of tti of the configurations of 
defects and will make a few comments about where our toy 
model goes wrong, compared to the full problem. In the next 
section, we will undertake a full and careful calculation of tti 
of the configurations of defects of a model with classifying 
space U(7V)/0(A^). 

A free fermion Hamitonian with no symmetry, but Umited 
to 8 bands, can be expanded about its minimum energy point 
in the Brillouin zone as 



H = ix{dili +niTi)x 



(42) 



From the considerations in the previous section, we learned 
that the space of mass terms is S"^, which is why we have 
written the mass term in the above form with a unit vec- 
tor n G S^. (More precisely, the space of mass terms is 
U(2)/0(2) = U(l) X 5'^, and the mass term can be written 
2^717273 j^^p^g-S7i7273 ^ Howcvcr, the extra U(l) plays no role 
here and can be ignored.) In the model of Section II, the three 
components of n correspond to dimerization in each of the 
three directions on the cubic lattice. In Teo and Kane's model 
they correspond to the real and imaginary parts of the super- 
conducting order parameter and the sign of the Dirac mass at 
a band inversion. 



TTrf-D-l (i?2+p-d) — TTo (Rp-D + l) 



(41) 



which is determined by the zero-th homotopy groups in Eq. 
(40). The important conclusion we obtain from this formula 
is that the classification of topological defect is determined by 
the dimension of the defect D and the symmetry class p, and 
is independent of the spatial dimension d. For p and D, we 
obtain the same 8x8 periodic table for the classification of 
topological defects. ^^ 

In the following we will focus on the 3D system with 
no symmetry, and discuss the generic cases in Sec. VII. 
Since a 3D system with no symmetry is classified by R-j = 
\J{N)/0{N), point-like defects in such a system are classified 
by 'K2{\i{N)/0{N)). In the stable limit, 7r2(U(A^)/0(7V)) = 
Z2. However, we note that, for smaller values of N be- 
low the stable limit, the classification is a little different, e.g. 
7r2(U(2)/0(2)) = 7r2(U(l) x S"^) = Z. The '8-band model' 
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FIG. 6: The pre-image of the North pole in the 5*^ target space of n is 
a collection of arcs connecting hedgehogs. A fixed tangent vector at 
the North pole defines a vector field along the arcs, thereby making 
them framed arcs, which may be viewed as ribbons. 
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FIG. 7: We will depict framed arcs or ribbons as arcs with twists 
accounted for by drawing kinks in the arcs, as shown above. 



Now consider defects in the n field, which are classified 
by TT2{S'^) = Z. Defects with winding number ±1 are pos- 
itive and negative hedgehogs. For simplicity, we will focus 
on these; higher winding number hedgehods can be built up 
from these. (In the real model, as opposed to the toy model, 
the 'hedgehogs' have a Z2 classification so there are no higher 
winding number hedgehogs and, in fact, they do not even have 
a sign.) As noted by Teo and Kane^*^, the n field around a 
hedgehog can be visualized in a simplified way, following 
Wilczek and Zee's discussion of the Hopf term in a 2 + 1- 
D 0(3) non-linear a model^^. The field n can be viewed as a 
map from the physical space where the electrons live, M^. If 
we assume that the total winding number is zero (equal num- 
bers of + and - hedgehogs) and that n approaches a constant 
at 00, we can compactify the physical space M^ so that it is 
S^. The target space of the map is 5^. The pre-image of the 
north pole A^ G 5^ is a set of arcs and loops. The choice of 
the north pole N is arbitrary, and any other point on the sphere 
would be just as good for the following discussion. Let's ig- 
nore the loops for the moment and focus on the arcs. Since 
n points in every direction at a hedgehog, the arcs terminate 
at hedgehogs. In fact, each arc connects a +1 hedgehog to 
a —1 hedgehog. We now pick an arbitrary unit vector in the 
tangent space of the sphere at N. This vector can be pulled 
back to S^ to define a vector field along the arcs which is 
clearly normal to the arcs. This is a framing, which allows us 
to define, for instance, a self-winding number for an arc. In- 
tuitively, we can think of a framing as a thickening of an arc 
into a ribbon. Thus, the field n allows us to to define a set of 
framed arcs connecting the hedgehogs - in other words, a set 
of ribbons connecting the hedgehogs. As the normal vector 
twists around an arc, the ribbon twists, as shown in Fig. 6 
(Although we will draw the ribbons as bands in the physical 
space, their width should not be taken seriously; they should 
really be viewed as arcs with a normal vector field.) 

Although these ribbons are strongly reminiscent of particle 
trajectories, it is important to keep in mind that they are not. 
A collection of ribbons connecting hedgehogs defines a state 
of the system at an instant of time. Ribbons, unlike particle 
trajectories, can cross. They can break and reconnect as the 
system evolves in time. As hedgehogs are moved, the ribbons 
move with them. 

A configuration of particles connected pairwise by ribbons 




FIG. 8: When two defects are exchanged, the n-field around them is 
modified. This is encapsulated by the dragging of the framed arcs as 
the defects are moved. 



is a seemingly crude approximation to the full texture defined 
by n. However, according to the Pontryagin-Thom construc- 
tion, as we describe in the next Section (and explain in Ap- 
pendix F), it is just as good as the full texture for topological 
purposes. Thus, we focus on the space of particles connected 
pairwise by ribbons. 

We now consider a collection of such particles and ribbons. 
For a topological discussion, all that we are interested in about 
the ribbons is how many times they twist, so we will not draw 
the framing vector but will, instead, be careful to put kinks 
into the arcs in order to keep track of twists in the ribbon, 
as depicted in Fig. 7. The fundamental group of their con- 
figuration space is the set of transformation which return the 
particles and ribbons to their initial configurations, with two 
such transformations identified if they can be continuously 
deformed into each other Consider an exchange of two +1 
hedgehogs, as depicted in Fig. 8. Although this brings the 
particles back to their initial positions (up to a permutation, 
which is equivalent to their initial configuration since the par- 
ticles are identical), it does not bring the ribbons back to their 
initial configuration. Therefore, we need to do a further mo- 
tion of the ribbons. By cutting and rejoining them as shown 
in Fig. 9a, a procedure which we call 'recoupling', we now 
have the ribbons connecting the same particles as in the initial 
configuration. But the ribbon on the left has a twist in it. So 
we rotate that particle by — 27r in order to undo the twist, as in 
Fig. 9b. 

Let us use ti to denote such a transformation, defined by 
the sequence in Figs. 8, 9a, and 9b. The f^s do not satisfy the 
multiplication rules of the permutation group. In particular, 
ti =/= t^^- The two transformations ti and t^^ are not dis- 
tinguished by whether the exchange is clockwise or counter- 
clockwise - this is immaterial since a clockwise exchange can 
be deformed into counter-clockwise one - but rather by which 
ribbon is left with a twist which must be undone by rotating 
one of the particles. 

To see that the operations ti, defined by the sequence in 
Figs. 8, 9a, and 9b, and t^^, defined by the sequence in 
10, are, in fact, inverses, it is useful to note that when they 
are performed sequentially, they involve two 27r twists of the 
same hedgehog. In 9b, it is the hedgehog on the left which is 
twisted; this hedgehog moves to the right in the first step of 10 
and is twisted again in the fourth step. One should then note 
that a double twist in a ribbon can be undone continuously by 
using the ribbon to "lasso" the defect, a famous fact related to 
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FIG. 9; (a) In order to restore the framed arcs so that they are con- 
necting the same defects, it is necessary to perform a recouphng by 
which they are reconnected. In order to keep track of the induced 
twist, it is easiest to perform the recoupling away from the overcross- 
ing. (b) The particle on the left must be rotated by —2n in order to 
undo a twist in the framed arc to which it is attached. 




FIG. 10: The sequence of moves which defines t^^ . (Here, the i* 
particle is at the top left and the (i + 1)"' is at the top right.) This 
may be contrasted with the sequence in Figs. 8, 9a, and 9b, which 
defines aj . 





FIG. 11: The sequence of moves which defines Xi: the defect on the 
left is rotated by 27r, the twist is transfered to the ribbon on the right 
by two recouplings, and then the defect on the right is rotated by 
-2n. (Here, the i"" defect is at the top left and the (i + 1)"' is at the 
top right.) The defects themselves are not moved in such a process. 



rotated by —2tt. Since a 47r rotation can be unwound, as de- 
picted in Fig. 16, xf = 1. 

Intuitively, one expects that Xi ~ ti^ since neither Xi nor ti^ 
permutes the particles and both of them involve 2tt rotations 
of both particles i and z + 1. To show that this is, in fact, 
the correct, we need to show that the history in Fig. 1 1 can 
be deformed into the sequence of Figs. 8, 9a, 9b repeated 
twice. If the history in Fig. 11 is viewed as a 'movie' and the 
sequence of Figs. 8, 9a, 9b repeated twice is viewed as another 
'movie' , then we need a one-parameter family of movies - or a 
'movie of movies' - which connects the two movies. We will 
give an example of such a 'movie of movies' shortly. With 
this example in hand, the reader can verify that Xi = ti^ by 
drawing the corresponding pictures, but we will not do so here 
since this discussion is superseded, in any case, by the the 
next section, where a similar result is shown for the A^ -^ 
oo problem by more general methods. We simply accept this 
identity for now. 

We now consider the commutation relation for the x^s. 
Clearly, for \i — i\ > 2, XiXj ~ XjXi. It is also intuitive 
to conclude that 



the existence of spin-1/2 and the fact that tti (5*0(3)) = Z2. 
This is depicted in Fig. 16 in Appendix E. It will be help- 
ful for our late discussion to keep in mind that ti not only 
permutes a pair of particles but also rotates one of them; any 
transformation built up by multiplying t^s will enact as many 
2tt twists as pairwise permutations modulo two. 

Thus far we have only discussed the +1 hedgehogs. We 
can perform the similar transformations which exchange —1 
hedgehogs. We will not repeat the above discussion for —1 
hedgehogs since the discussion would be so similar; further- 
more, in the iV — > 00 model which is our main interest, de- 
fects do not carry a sign, so they can all be permuted with each 
other. 

We have concluded that ti ^ t^^ and, therefore, the group 
of transformations which bring the hedgehogs and ribbons 
back to their initial configuration is not the permutation group. 
This leaves open the question: what is tj^? The answer is 
that tf can be continuously deformed into a transformation 
which doesn't involve moving any of the particles - Teo and 
Kane's 'braidless operations'. Consider the transformation Xi 
depicted in Fig. 11. Defect i is rotated by 27r, the twist is 
transferred from one ribbon to the other, and defect z + 1 is 



•^i-^i-\-l •^i-\-l'^i 



(43) 



since the order in which twists are transferred is seemingly 
unimportant. However, since this is a crucial point, we ver- 
ify it by showing in Figure 12 that the sequence of moves 
which defines XiXi+i (a 'movie') can be continuously de- 
formed into the sequence of moves which defines Xi^iXi (an- 
other 'movie'). Such a deformation is a 'movie of movies'; 
going from left-to-right in Fig. 12 corresponds to going for- 
ward in time while going from up to down corresponds to de- 
forming from one movie to another. 

Thus, we see that the equivalence class of motions of the 
defects (i.e. tti of their configuration space) has an Abelian 
subgroup generated by the .t^s. Since xf = 1 and they all 
commute with each other, this is simply n — 1 copies of Z2, 
or, simply, (Z2)"~^. 

In order to fully determine the group of transformations 
which bring the hedgehogs and ribbons back to their initial 
configuration, we need to check that the i^s generate the full 
set of such transformations - i.e. that the transformations de- 
scribed above and those obtained by combining them exhaust 
the full set. In order to do this, we need the commutation 
relations of the t^s with each other. Clearly, titj = tjti for 
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FIG. 12: The sequence of moves which defines XiXi-^-i is shown in the top row. The sequence of moves which defines 0:^+13::^ is shown in 
the bottom row. The rows in between show how the top row can be continuously deformed into the bottom one. Such a deformation of two 
different sequences is a 'movie of movies' or a two-parameter family of configurations. Moving to the right increases the time parameter while 
moving down increases the deformation parameter which interpolates between xtXi+i and Xi+\Xi. 



\i — j\ > 2 since distant operations which do not involve the 
same hedgehogs nor the same ribbons must commute. On the 
other hand operations involving the same hedgehogs or rib- 
bons might not commute. For instance, 



^«*^i+l ■^i-^i'\-lf'i 



(44) 



To see why this is true, note that if we perform x^+i first, then 
defects i + 1 and i + 2 are twisted by 2tt. However, ti then 
permutes i and i + 1 and twists i by 27r. Thus, the left-hand- 
side permutes i and i + 1 and only twists i + 2. The (i + 1)**^ 
hedgehog was twisted by Xi^i and then permuted by ti so that 
it ended up in the i*"^ position, where it was twisted again in 
the last step in t^; two twists can be continuously deformed to 
zero, so this hedgehog is not twisted at all. The right-hand- 
side similarly permutes i and i + 1 and only twists i + 2 by 
2tt. The reader may find it instructive to flesh out the above 
reasoning by constructing a movie of movies. 

The multiplication rule which we have just described (but 
not fully justified) is that of a semi-direct product, which is 
completely natural in this context: when followed by a per- 
mutation, a transfer of twists ends up acting on the permuted 
defects. The twists Xi form the group (Z2)"~^ which we can 
represent by n-component vectors all of whose entries are 
or 1 which satisfy the constraint that the sum of the entries is 
even. The entries tell us whether a given hedgehog is twisted 
by 27r or not. In any product of x^s, an even number of hedge- 
hogs is twisted by 27r. Now consider, for n odd, the group 
elements given by 

"■j — Xn-lXn-3 ■ ■ ■ XiJ^3Xi^iXi^2Xi^4 . . . Xiti (45) 

for i odd and 



X'n—lX'n—S • • • XiJ^2XiXi^\Xi 



■XlU 



(46) 



for i even. From (44), we see that ai^ = 1. The group element 
(Ti permutes the i^^ and (i + 1)*'^ hedgehogs and twists all of 
the hedgehogs. Thus, the aiS generate a copy of the permuta- 
tion group Sn- The cr^s do not commute with the x^s, however; 
instead they act according to the semi -direct product structure 



noted above. On the other hand, the situation is a bit different 
for n even. This may be a surprise since one might expect that 
n even is the same as n odd but with the last hedgehog held 
fixed far away. While this is true, exchanging the last hedge- 
hog with the others brings in an additional layer of complexity 
which is not present for n odd. The construction above, Eqs. 
45, 46, does not work. One of the hedgehogs will be left un- 
twisted by such a construction; since subsequent aiS will per- 
mute this untwisted hedgehog with others, we must keep track 
of the untwisted hedgehog and, therefore, the (TiS will not gen- 
erate the permutation group. In the even hedgehog number 
case, the group of tranformations has a (^2)"^^ subgroup, as 
in the odd case, but there isn't an Sn subgroup, unlike in the 
odd case. To understand the even case, it is useful to note that 
in both cases, every transformation either (a) twists an even 
number of ribbons, which is the subgroup (Z2)"~^; (b) per- 
forms an even permutation, which is the subgroup A„ of Sn', 
or (c) twists an odd number of ribbons and performs an odd 
permutation. Another way of saying this is that the group of 
transformations is the 'even part' of (Z2)" x Sn'- the subgroup 
of (Z2)" XI Sn consisting of those elements whose (Z2)" par- 
ity added to their 5„ parity is even. In the odd hedgehog num- 
ber case, this is the semidirect product (22)"^^ x Sn', in the 
even hedgehog number case, it is not. As we will see in Sec- 
tion VI, the difference between the even and odd hedgehog 
number cases is related to the fact that, for an even number of 
hedgehogs, the Hilbert space decomposes into even and odd 
total fermion number parity sectors. By contrast, the situation 
is simpler for an odd number of hedgehogs, where the parity 
of the total fermion number is not well-defined and the repre- 
sentation is irreducible. 

To summarize, we have given some plausible heuristic ar- 
guments that the 'statistics' of +1 hedgehogs in a model 
of 2N — 8 Majorana fermions is governed by a group 
£{{^2)"- XI Sn), the 'even part' of (Z2)" xi 5„: those elements 
of (Z2)" XI Sn in which the parity of the sum of the entries of 
the element in (Z2)" added to the parity of the permutation 
in Sn is even. (The same group governs the —1 hedgehogs). 
Rather than devoting more time here to precisely determining 
the group for the toy model, we will move on to the prob- 
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lem which is our main concern here, a system of 2N -^ oo 
Majorana fermions. This problem is similar, with some im- 
portant differences. (1) The target space is no longer 5*^ but 
is, instead, U{N)/0{N). (2) Consequently, the defects do 
not carry a sign. There is no preferred pairing into ± pairs; 
the defects are all on equal footing. All 2n of them can be 
exchanged. (3) The group obtained by computing tti of the 
space of configurations of 2n defects then becomes the direct 
product of the 'ribbon permutation group' T^i with a trivial Z, 
72n = Z X T^n- The ribbon permutation group 75„ is given 
by Ta^ = Z2 X S((Z2)2« X 52„), where ^((Za)^" x 52„) is 
the 'even part' of (Z2)^" x S'2„. 



V. FUNDAMENTAL GROUP OF THE MULTI-DEFECT 
CONFIGURATION SPACE 

In Section III we concluded that the effective target space 
for the order parameter of a system of fermions in 3D with no 
symmetries is U{N)/0{N) - which, as is conventional, we 
will simply call U /O, dropping the N in the large-A^ limit. 
This enables us to rigorously define the space of topological 
configurations, K2n, of 2n hedgehogs in a ball, and calculate 
its fundamental group TTi{K2n), thereby elucidating Teo and 
Kane's^*^ hedgehog motions and unitary transformations. 

We now outline the steps involved in this calculation: 

• We approximate the space U/O by a cell complex 
(or CW complex), C, a topological space constructed 
by taking the union of disks of different dimensions 
and specifying how the boundary of each higher- 
dimensional disk is identified with a subset of the lower- 
dimensional disks. This is a rather crude approxima- 
tion in some respects, but it is sufficient for a homotopy 
computation. 

• We divide the problem into (a) the motion of the hedge- 
hogs and (b) the resulting deformation of the field 
configuration between the hedgehogs. This is accom- 
plished by expressing the configuration space in the fol- 
lowing way. Let us call the configuration space of 2n 
distinct points in three dimensions X2n- (For the sake 
of mathematical convenience, we will take our physi- 
cal system to be a ball B"^ and stipulate that the points 
must lie inside a ball B^. Let's denote the space of 
field configurations by A^2n- This space is the space 
of maps to U/O from B^ with 2n points (at some stan- 
dard locations) excised. The latter space is denoted by 
B^ \ 2n standard points. Since we will be approximat- 
ing U/O by C, we can take A^2n to be the space of 
maps from B^ \ 2n standard points to C with boundary 
conditions at the 2n points specified below. Then, there 
is afibration of spaces: 



M2n^K2n 
X2n 



• We introduce another two fibrations which further di- 
vide the problem into more manageable pieces: 



R2n -^ K2n 

Y2n 



J^2n^Y2n 
X2n 

The original fibration is kind of a "fiber-product" of the 
two new fibrations. Here, i?2n is essentially the space 
of order parameter textures interpolating between the 
hedgehogs, and y2n is the space of configurations of 2n 
points with infinitesimal spheres surrounding each point 
and maps from each of these spheres to C. N2n is the 
space of maps from 2n infinitesimal spheres to C, with 
each one of the spheres surrounding a different one of 
the 2n points (at some standard locations) excised from 
B^. We call these order parameter maps from infinites- 
imal spheres to C "germs". 

• Having broken the problem down into smaller pieces 
by introducing these fibrations, we use the fact that a 
fibration F ^ E -^ B induces a long exact sequence 
for homotopy groups 

. . . ^ n,{E) -^ TT,{B) -^ 7r,_i(F) ^ TT,-i{E) -^ ... 

For instance, applying this to the fibration Ai2n ^ 
K2n -^ X2n Icads to the exact sequence ... -^ 

7ri(Al2n) -^ '^l{K2n) ^ T^l{X2n) ^ 1- It folloWS 

that TTi{K2n) is an extension of the permutation group 
S2n ~ T^i{X2n)- By itsclf, the above long exact se- 
quence is not very helpful for computing any of the 
homotopy groups involved unless we can show by in- 
dependent means that two of the homotopy groups are 
trivial. Then the homotopy groups which lie between 
the trivial ones in the sequence are tightly constrained. 

• We directly compute that 7ri(A/'2n) = (Z2)^" and 
T^i{X2n) — S2n- We show that the homotopy exact 
sequence then implies that 'Ki{Y2n) = (Z2)^" x 52„. 

• We compute the homotopy groups of i?2n, defined by 
the fibration i?2n — > ^2n — > Y2n- This computation 
involves a different way from the cell structure of think- 
ing about the topology of a space, called the "Postnikov 
tower", explained in detail Appendix E. The basic idea 
is to approximate a space with spaces with only a few 
non-trivial homotopy groups. (This is analogous to the 
cell structure, which has only a few non-trivial homol- 
ogy groups.) The simplest examples of such spaces are 
Eilenberg-Mac Lane spaces, which only have a single 



20 



non-trivial homotopy group. The Eilenberg-Mac Lane 
space K{A, m) is defined for a group A and integer m 
as the space with homotopy group Trm{K{A, m)) = A 
and TTkiKiA, m)) = for all k ^ m. (The group A 
must be Abelian for to > 1.) Such a space exists and is 
unique up to homotopy. A space T with only two non- 
trivial homotopy groups can be constructed through the 
fibration K{B, n) ^ T ^ K{A, m). The space T has 
T^m{T) ~ A and 7r„(T) = B, as may be seen from 
the corresponding long exact sequence for homotopy 
groups. Continuing in this fashion, one can construct 
a sequence of such approximations M„ to a space M. 
They are defined by 7rfc(M„) — TTk{M) for k < n and 
T^k{Mn) — for k > n. They can be constructed it- 
eratively from the fibration K{A, n) -^ M„ -^ A'/„_i, 
where 7r„(M) = A. 

• With T:i{Y2n), 71-2 (i^2«), 7i"o(i?2«) and 7ri(i?2„) in hand, 
we compute the desired group tti {K2n) from the homo- 
topy exact sequence. 

We now go through these steps in detail. 

Approximating U/O by a cell complex. Depending on 
microscopic details, gradients in the overall phase of the 
fermions may be so costly that we wish to consider only con- 
figurations in which this overall phase is fixed. We will re- 
fer to this as the scenario in which 'phase symmetry is bro- 
ken'. In this case, the effective target space is SU /SO, the 
non-phase factor of U/O ^ U{1)/0{1) x SU/SO. In this 

case, we simplify matters by replacing SU/SO by U/O, the 

universal cover of U/O. U/O is homotopy equivalent to 
SU/SO, so this substitution is harmless. This substitution 
results in a reduced configuration space K2n and we will 
concentrate first on calculating 'ni{K2n)- In an appendix, 
we show that this reduction essentially makes no difference: 

^l(i^2n)=7ri(X2«)xZ. ^_^ 

We now define a cell complex C approximating U/O. 
In constructing this cell structure, we are not interested in 

the beautiful homogeneous nature of U/O but rather only 
its homotopy type. The homotopy type of a space tells 
you everything you will need to know to study deformation 
classes of maps either into or out of that space. An im- 
portant feature of any homotopy type is the list of homo- 
topy groups (but these are by no means a complete charac- 
terization in general). For U/O, the homotopy groups are 
T^iiujd) = 0,Z2,Z2,0,Z,0,0,0,Z for i = l,...,9and 

thereafter TTi{U /O) cycles through the last eight groups. (For 
U/O, the first group would be Z.) 

Because U/O is simply-connected, but has nontrivial 7r2, it 
natural in building a cellular model for its homotopy type to 
begin with S"^ . Since tt2{S'^) — Z and we only need a Z2 for 

TT2(U /O), we should kill off the even elements by attaching 
a 3-cell D^ using a degree-2 map of its boundary 2-sphere 
to the original S^. For future reference, take this map to be 
[0, (p) -^ (26, (p) in a polar coordinate system where the north 



pole N — (tTjO). Similarly, a 4-cell is attached to achieve 

TT^{U /O) = Z2. The necessity of the 4-cell is proved (Fact 1) 
below. 

The preceding logic leads us to the cell structure: 

C^S^ IJ D^ [j £!'' IJ cells of dimension > 5 (47) 

degree— 2 2Hopf 

Since we are only trying to compute the fundamental group 
TTi{K2n) from our various homotopy long exact sequences, 
we do not have to figure out the higher cells (dimension > 5) 

of U/O. We will, however, verify that tt^{U/0) is generated 

by the Hopf map into the base 5*^ C U/O. 

To summarize, we will henceforth assume that the order 
parameter takes values in the cell complex C. Although C 
is a crude approximation for U/O, it is good enough for the 
topological calculations which follow. 

Dividing the problem into the motion of the hedgehog 
centers and the deformation of the field configuration. Let 
us assume that our physical system is a ball of material B^. 
Let n > be the number of hedgehog pairs in the system. 
A configuration in K2n is a texture in the order parameter, 
^{x) : B^ ^ C, which satisfies the following boundary con- 
ditions at the boundary of B^ and at the 2n hedgehog loca- 
tions (which are singularities in the order parameter). The 
order parameter has winding number at the boundary of the 
ball, dB^ and winding number 1 around each of the hedgehog 
centers. (Recall that tt2{C) = Z2, so the winding number can 
only be or 1). 

From its definition, K2n is the total space of a fibration: 

M2n^K2n 
X2n 

The above diagram suggests that we should think of the fi- 
bration Ai2n ^ K2n -^ ^2n in the following way: above 
each point in X2n there is a fiber Al2n; the total space formed 
thereby is K2n- (This is not quite a fiber bundle, since we 
do not require that there be local coordinate charts in which 
K2n is simply the direct product.) Here, X2n is the simply 
the configuration space of 2n distinct points in B^. We write 
this formally as X2n = Yiili ^^ \ big diagonal. (The big di- 
agonal consists of 2n-tuples of points in B^ where at least two 
entries are identical.) The space A^2n consists of maps from 
B^ \ 2n points in a fixed standard position to C with the pre- 
scribed winding numbers given in the preceding paragraph. 

Germs of order parameter textures. It is helpful to in- 
troduce an intermediate step in the fibration. Define a point in 
Y2n as a configuration in X2n together with a "germ" of $(x), 
which we call $(x), defined only near dB'^ and the 2n points. 
The idea behind the germ $(x) is to forget about the order pa- 
rameter $(a:) except for its behavior in an infinitesimal neigh- 
borhood around each hedgehog center and at the boundary of 
the system. $(a;) must satisfy the same boundary conditions 
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FIG. 13: B"^ \2n points in standard position. The space A42n con- 
sists of maps from this manifold to C. 



as ^(x) itself. We take $(a;) to be constant on dB^ and to 
have winding number 1 around each of the hedgehog centers. 
With this definition, we now have the fibration: 



I 
X2n 

where M2n is the space of (germs of) order parameter textures 
<I> from the neighborhoods of the 2n fixed standard points and 
dB^ to C. We will henceforth replace discussion of germs 
with the equivalent and simpler concept of maps on dB^ U 



(u^ 



standard point. Thus, 



Sf ) where S'^ is a small sphere surrounding the zth 



(48) 



J\f2n C Maps (as^ \j[j SJ\^C 



We now define Q2n as the ball B^ with small balls (denoted 
below by interior(5|)) centered about the hedgehogs deleted: 



Q2 



B^\\J interior(S'f 



(49) 



i=l 



for fixed standard positions i = 1, • • • , 2n. Then i?2n is the 
space of order parameter textures on Q2n which satisfy the 
boundary condition that the winding number is on dB^ and 
1 on each of the small spheres. With this definition, we have 
the fibration: 



R2n -^ K2, 



Y2, 



Given the cell structure C, we can specify the boundary con- 
ditions for the order parameter precisely. On dB^, the or- 
der parameter is equal to the North Pole in 5*^. (Recall that 
S"^ C C is the bottom cell of the structure C which we are us- 
ing to approximate U/O.) On each of the spheres Sf, the order 
parameter $(a;) defines a map from Sf -^ S"^ which is the 
identity map (where, again 5*^ is understood as a subset of the 



order parameter space S'^ C C). This ensures that the order 
parameter has the correct winding numbers at the boundaries 
of Q2n- In essence, what we have done in writing Eq. V is 
to break up an order parameter texture containing hedgehogs 
into (a) the hedgehogs together with the order parameter on in- 
finitesimal neighborhoods around them (i.e. 'germs') and (b) 
order parameter textures in the intervening regions between 
the hedgehogs. The space of configurations (a) is Y2n', the 
space of configurations (b) is i?2n- 

The name i?2n is for "ribbons." As we saw in Section IV, if 
the order parameter manifold were 5^, we could summarize 
an order parameter texture by looking at the inverse image of 
the North Pole iV C 5^ and a fixed tangent vector at the North 
Pole. The inverse images form a collection of ribbons. Now, 
the order parameter manifold is actually (approximated by) C, 
but the bottom cell in C is S"^. The effect of the 3 -cell is that 
hedgehogs lose their sign, so there is no well-defined "arrow" 
running lengthwise along the ribbons. The 4-cell allows the 
"twist" or framings of ribbons to be altered at will by ±2. 

Long exact sequence for homotopy groups. It is very 
convenient to use fibrations to calculate homotopy groups. 
(For those interested in K-Xhtory, the last two chapters of 
Milnor's Morse Theory^^ are a must read and exhibit these 
methods with clarity.) As noted above, fibrations have all the 
homotopy properties of fiber bundles but are (often) found 
arising between function spaces where it would be a lot of 
work - and probably a distraction from important business - 
to attempt to verify the existence of locally trivial coordinate 
charts. Operationally, fibrations share with fiber bundles the 
all-important "homotopy long exact sequence": 



F 



E 

B 

we have: 

> 7r,+i(B) ^ ^,{F) ^ t:,{E) ^ 7r,(B) ^ 7r,_i(i^) ^ • 

We now compute Tii{Y2n) from the exact sequence: 
d d 

■^2 (^2« ) > TTi {M2n ) -^ TTi {J2n) -^ T^\ (^2n) > TTQ (A2n) 

(50) 
We can compute two of the homotopy groups in this equa- 
tion by inspection, tti (X2„) is clearly the symmetric group of 
point exchange: 



7ri(X2 



^2 



(51) 



Meanwhile, 7ri(A/2n) amounts to (products of) loops of maps 
from the Sf to C and reduces to 2n copies of the third ho- 
motopy group of C (and, therefore, to ■k-3,{U jO)). Thus, 

T^l{M2n) - (^2)2": 



2n 



7ri(AA2„)-[]^i(Maps(5f,(7/0)) = J] ^3(C//0) 

i=:l In copies 
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(Z2) 



2n 



(52) 



To proceed further, we need to evaluate boundary maps in 
the homotopy exact sequence. In Appendix D, we explain 
boundary maps through the example of the Hopf map. Con- 
sider Eq. 50. 7r2(X2„) is generated by the ( ^) different 2- 
parameter motions in which a pair of hedgehogs come close 
together and explore the 2-sphere of possible relative posi- 
tions around their center of mass. This 2-parameter family of 
motions involves no "rotation" of the maps $ which describe 
7ri(A/2n) (i-e the order parameter configuration in the neigh- 
borhood of each hedgehog does not rotate as the hedgehogs 
are moved), so the left most d map in Eq. 50 is zero. Sim- 
ilarly, a simple exchange of hedgehogs produces no twist of 
the order parameter configuration in the neighborhood of ei- 
ther hedgehog, so the second d map of Eq. 50 is also zero. 
Thus, we have a short exact sequence: 



1 



^r- 



■MY2n)- 



(3 



>^2 



1 



To derive this short exact sequence, we used the triviality of 
the boundary maps noted above and Eqs. 51, 52 to simplify 
Eq. 50. 

There is a natural group homomorphism s: 

S : S2n^ 7ri(F2«) 

which associates to each permutation a motion of hedgehogs 
which permutes the hedgehogs in Y2n but does not rotate the 
order parameter configurations $ near the hedgehogs. Then 
j3 o s = zdsjn ■ ^ri other words, the sequence is split: 



First consider di. The kernel of di is represented by loops 
in Y2n which extend to loops in i?2n. A loop 7 in l2n is a 
motion of the hedgehogs together with rotations (about the 
spheres Sf) of $ which brings the system back to its initial 
configuration. If a loop 7 is in kernel of di, then there is a 
corresponding loop in the configuration space of ribbons in 
B^ (obtained by lifting 7 to K2n)- 

The next steps are to compute ttq {R2n ) and tti (i?2ri ) ■ These 
computations are detailed in Appendix E, where we see that 

using the cell structure C which we introduced for U /O is 
tricky as a result of the higher cells. Thus, we instead intro- 
duce the "Postnikov tower" for U /O which allows us to make 
all the calculations we need. We find that 7ro(i?2n) = ^2 and 

^l(i?2n) - (Z2)2". 

Thus, T^i{K2n) sits in the following exact sequence: 
d ~ d 

7r2 iY2n) >■ TTl (i?2n) ^ TTi {K2„) ^ T^l (^2n) ^ TTo (i?2n) 



(Z2)2" 



^11 l\\ 

(Z2)2« X S2n Z2 



Recall that when we studied TT2{Y2n), we found (exactly as 
in the case of TT2{X2n)) that there are (^^) generators corre- 
sponding to relative 2-parameter motions of any pair of hedge- 
hogs around their center of mass. This can be used to un- 
derstand the map 82 : 7r2(l2n) ^^ T^2{R2n)- The image of 
any center of mass 2-motion is the "bag" containing the corre 
sponding pair of hedgehogs. Thus, coker(92) — ^2', it is 
modulo the even sublattice (vectors whose coordinate sum is 
zero in Z2). Thus, we have the short exact sequence: 



2n 
2 



/3 



■Zi" 



■7ri(y2n): 



-> 1 



1 )-coker((92) 

Z2 



-7ri(if2„)^ker(ai) 



1 



(54) 



Thus, TTi (Y2n) is a semi-direct product. To determine tti (^2™) 
completely, it only remains to identify how s{S2n) acts on the 
twist factors Zg" under conjugation. It is quite clear that this 
action is the only natural one available: s{p) acts on Z2" by 
applying the permutation p to the 2n coordinates of Z|". So, 



Tri{Y2; 



Zi" 



>! S2n with group law: 



{v,p)o {v',p') = {v+p{v'),pop') 



(53) 



where v € Z2" is a Z2 -vector, p e S2n a permutation, and 
p{v') the natural action of S'2„ on Z2" applied to v' . Note 
that this is precisely the multiplication rule which we obtained 
pictorially in Section IV. 

Computing the homotopy groups of i?2n, the space 
of order parameter textures interpolating between the 
hedgehogs. Of course, computing 'ni{Y2n) only gets us part 
of the way home. Our ultimate goal is to compute TTi{K2n)- 
Thus, we now turn to the homotopy long exact sequence: 



T^2{Y2, 



-^ 7ri(i?2„) -^ 7ri(ii'2„) -> TTl{Y2n) — ^ 7ro(i?2«) 



The kernel ker(9i ) consists of the part of tti (Y) associated 
with even (2tt) twisting. As shown in Figure 9, a simple ex- 
change is associated to a total twisting of ribbons by ±27r. 

Thus, di{v,p) = X]j=i ^i + parity(p) € Z2. We use the 
notation E{Z^ xi S,n) forker(ai). 

Note: If 771 = 2, Z™ xi Sm is the dihedral group D4 and 
its "even" subgroup ker(9i) ^ Z4. This shows that for m 
even, the induced short exact sequence does not split, and the 
extension is more complicated: 



1 



zi 



E{Zf X ^2,0 



S2 



1 



There is a final step required to solve the extension problem 
54 and finish the calculation of TTi{K2n)- We geometrically 
construct a homomorphism s : ker(9i) = E{1j'^'' x S'2„) -^ 
TTi{K2n) which is a left inverse to the projection. 

This will show that TTi{K2n) is a semidirect product Z2 x 
ker(9i), but since Z2 has no nontrivial automorphism, the 
semidirect product is actually direct: 
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^l(X2n)=Z2XS(Z2">^^2«) 



(55) 



To construct s, note that all elements of ker(9i ) can be re- 
alized by a loop 7 of maps into the bottom 2-cell of Eq. 47 
5^. Still confining the order parameter (map) to lie in S*^, 
such a loop lifts to an arc 7 of ribbons representing an arc 
in K2n- We may choose the lift so that as the ribbons move, 
they never "pass behind" the 2n hedgehogs. (For example, we 
may place the hedgehogs on the sphere of radius = ^ inside 
the 3-ball B"^ (assumed to have radius = 1) and then keep 

all ribbons inside B\ . These arcs may be surgered (still as 

2 

preimages of A^ C 5*^) so that they return to their original po- 
sition except for a possible accumulation of normal twisting 
t2Tr. Since 7 G ker(9i), t must be even. Now, allowing the 
order parameter (map) to leave S*^ and pass over the 4-cell (of 
Eq. 47), attached by 2Hopf : S^ — >■ S'^, we may remove these 
even twists. (The 4-cell can introduce small closed ribbons 
with self-linking = 2 in a small ball. These small ribbons can 
be surgered into other ribbons.) This lifts a generating set of 
ker(9i) into TTi{K2n) as a set theoretic cross section (left in- 
verse to tt). But what about relations? Because the entire loop 

is constant outside B\ , the corresponding homology class in 

2 

i?2(Q X /;Z2) = 7ri(i?) is trivial, so s is actually a group 
homomorphism. 



VI. REPRESENTATION THEORY OF THE RIBBON 
PERMUTATION GROUP 

In this section, we discuss the mathematics of the group %n 
and its representation. The purpose of this section is to show 
that a direct factor of Tm, called the even ribbon permutation 
group, is a ghostly recollection of the braid group and the Teo- 
Kane representation of the even ribbon permutation group is a 
projectivized version of the Jones representation of the braid 
group at a 4*'^ root of unity, i.e. the representation relevant to 
Ising anyons. 



A. Teo-Kane fundamental groups 

In Section V, we consider only even number of hedgehogs 
for physical reasons. In this section, we will include the odd 
case for mathematical completeness. 

The Teo-Kane fundamental group is the fundamental group 
of the Teo-Kane configuration space K^. As computed in 
Section V, T™ == T:i{K.m) = Z x Z2 x E{1'^ x Sm), where 
the subgroup T^^ = Z2 x ^^(Z™ x Sm) is called the ribbon 
permutation group. Here, E{Z^ x S„i) is the subgroup of 
Z™ X Sm comprised of elements whose total parity in Z™ 
added to their parity in Sm is even. In the following, we will 
call the group Gm = E{Z'!^ x Sm) the even ribbon permu- 
tation group because it consists of the part of tti (Ym) associ- 
ated with even (27r) twisting. For the representations of the 
Teo-Kane fundamental groups, we will focus on the even rib- 



bon permutation groups Gm- No generality is lost if we con- 
sider only irreducible representations projectively because ir- 
reducibles of Z and Z2 contribute only overall phases. But 
for reducible representations, the relative phases from repre- 
sentations of Z and Z2 might have physical consequences in 
interferometer experiments. 

The even ribbon permutation group Gm is an index= 2 sub- 
group of Z™ X Sm ■ To have a better understanding of Gm, we 
recall some facts about the important group Z™ x Sm- The 
group Z™ X Sm is the symmetry group of the hypercube Z™, 
therefore it is called the hyperoctahedral group, denoted as 
Gm- Cm is also a Coxeter group of type Bm or Gm, so in 
the mathematical literature it is also denoted as Bm or BCm- 
To avoid confusion with the braid group Bm, we choose to 
use the hyperoctahedral group notation Gm- The group Gm 
has a faithful representation as signed permutation matrices 
in the orthogonal group 0{m): matrices with exactly one 
non-zero entry ±1 in each row and column. Therefore, it 
can also be realized as a subgroup of the permutation group 
52m, called signed permutations: a : {±1, ±2, • • • , ±m] -^ 
{±1, ±2, • • • , ±.m] such that cr(— i) ~ — cr(i). 

We will denote elements in Gm by a pair (6,5), where 
h = (bi) G Z™ and g G Sm- Recall the multiplication of 
two elements (b, g) and (c, h) is given by (6, g) ■ (c, h) — 
{b + g.c,gh), where g.c is the action of g on c by per- 
muting its coordinates. Let {e^} be the standard basis el- 
ements of M™. To save notation, we will also use it for 
the basis element of Z™. As a signed permutation matrix 
Ci introduces a —1 into the i*'* coordinate Xi. Let Si be 
the transposition of Sm that interchanges i and i + 1. As 
a signed permutation matrix, it interchanges the coordinates 
Xi,Xi+i. There is a total parity map det : Gm —> ^2 de- 
fined as det{b,g) — '^"Li^i + parity (g) mod 2. We de- 
note the total parity map as det because in the realization of 
Gm as signed permutation matrices in 0(m), the total par- 
ity is just the determinant. Hence Gm, as the kernel of det, 
can be identified as a subgroup of SO{m). The set of ele- 
ments ti — (ei, Si), i = 1, • • • ,m — 1 generates Gm- As a 
signed permutation matrix, ti{xi, ■ ■ ■ , Xi,Xi+i, ■ ■ ■ , Xm) ~ 

Given an element {b,g) G Gm, let Z™^^ be identified as 
the subgroup of Z™ such that X]t=i ^i is even. Then we have: 

Proposition VI.l. 1. For m > 2, the even ribbon per- 
mutation group Gm has a presentation as an abstract 
group 

<ti,--- ,tm-l\tf = l,i = !,•• • ,TO- 1, 

(W+i)^ = (ti'U+i)^ = l,i = 1, • • • ,TO - 2 > . 

2. The exact sequence 

1 ^ ^m-l ^ /^ ^ e s 1 

1 ? ili2 7 ^m ^ "-"m ^ -•- 

splits if and only ifm is odd. 

3. When m is even, a normalized 2-cocycle f(g, h) : Sm x 



o„ 



ZI; 



for the extension of Sm by Z™ above 



24 



can be chosen as f{g,h) — if g or h is even and 
f{g, h) = 65,(1) if g and h are both odd. 

We briefly give the idea of the proof of Prop. VI. 1. For (1), 
first we use a presentation of Cm as a Coxeter group of type 
Bra- < ri,--- ,r„i|rf = 1, (rir2)^ = 1, (rjri+i)^ = 1,« = 
2, • • • ,TO — 1 >. Then the Reidemeister-Schreier method^"* 
aUows us to deduce the presentation for Gm above. For (2), 
when m is odd, a section s for the spHtting can be defined 
as s{g) = (0,5) if g is even and 3{g) = ((11 • • • l),.g) if 
g is odd. When m is even, that the sequence does not split 
follows from the argument in^^. For (3), we choose a set map 
s{g) — (0, g) if g is even and s{g) — (ei , g) if g is odd. Then 
a direct computation of the associated factor set as on page 91 
of^*" gives rise to our 2-cocycle. 

As a remark, we note that there are another two obvious 
maps from the hyperoctahedral group C„i to Z2. One of them 
is the sum of bits in b of (6, g). The kernel of this map is the 
Coxeter group of type Dm, which is a semi-direct product of 
Z™^^ with Sm- The two groups Dm and Gm have the same 
order, and are isomorphic when m is odd (the two splittings 
induce the same action of Sm on Z™~^), but different when 
m is even. To see the difference, consider the order 2 auto- 
morphism (f) : Gm -^ Cm given by (t>{x) = det{x)x. Its 
restriction is the identity on Cm, but non-trivial on Dm- 



B. Teo-Kane unitary transformations 



Suppose there are to hedgehogs in i? . A unitary transfor- 
mation Ti 



«j 



nii'ij of the Majorana fermions is associated 
with the interchange of the (i,j)-pair of the hedgehogs. Inter- 
changing the same pair twice results in the "braidless" opera- 
tion T? ~ lilj- The Majorana fermions 7^,? — 1,2,--- ,?n 
form the Clifford algebra CImi'C). Therefore, the Teo-Kane 
unitary transformations act as automorphisms of the Clifford 
algebra 



T 



7 



TtjjT'^ 



The projective nature of the Teo-Kane representation rears its 
head here already as an overall phase cannot be constrained by 
such actions. Do these unitary transformations afford a lin- 
ear representation of the Teo-Kane fundamental group Km"^ 
If so, what are their images? The surprising fact is that the 
Teo-Kane unitary transformations cannot give rise to a linear 
representation of the Teo-Kane fundamental group. The re- 
sulting representation is intrinsically projective. We consider 
only the even ribbon permutation group Gm from now on. 

To define the Teo-Kane representation of Cm, we use the 
presentation of Gm hy ti,i — 1,- ■■ ,m — I in Prop. VI. 1. 
The associated unitary matrix for ti comes from the Teo- 
Kane unitary transformation T^ i+i. As was alluded above, 
there is a phase ambiguity of the Teo-Kane unitary matrix. 
We will discuss this ambiguity more in the next subsection. 
For the discussion below, we choose any matrix realization 
of the Teo-Kane unitary transformation with respect to a ba- 
sis of the Clifford algebra GlmiC). A simple computation 
using the presentation of Gm in Prop. VI. 1 verifies that 



the assignment of Ti.i+i to ti indeed leads to a represen- 
tation of Gm.- Another verification follows from a relation 
to the Jones representation in the next subsection. We can 
also check directly that this is indeed the right assignment for 
T^i+i ■■ 7i -^ -7i,7i+i -^ -li+i- In the CUfford alge- 
bra Glm (C), ji , 7i+i correspond to the standard basis element 



■-i 1 ^i-\- 



1. The element t^ of Gm is (e^ 



"t+ 



1, 1). As a signed 



permutation matrix, tf sends {xi, ■ ■ ■ ,Xi, x^+i, • • • , Xm) to 
{xi, ■ ■ ■ , —Xi, — Xi+i, • • • , Xm), which agrees with the action 
ofTi%iOn7,,7j+i. 

To see the projectivity, the interchange of the {i,i + l)-pair 
hedgehogs corresponds to the element ti = (ei,Si) G Gm- 
Performing the interchange twice gives rise to the element 
tf = {ci + Gi+i, 1), denoted as Xi. Since x/s are elements 
of a subgroup of Gm isomorphic Z™^^, obviously we have 
XiXi+i — Xi+iXi. On the other hand, T^-^i+iTj^^ ^^3 = 
Jzli+i ■ 7j+i7*+2 = 7j7»+2, and T^+i^^+2Tli+i = lt+ili+2 ■ 
7i7i+i = ~7i7i+2- Since it is impossible to encode the —1 in 
the Xi's of Cm, the representation has to be projective. Note 
that an overall phase in Tij will not affect the conclusion. In 
the next subsection, we will see this projective representation 
comes from a linear representation of the braid group — the 
Jones representation at a 4*^* root of unity and the — 1 is en- 
coded in the Jones-Wenzl projector p3 = 0. 

To understand the images of the Teo-Kane representation 
of Cm, we observe that the Teo-Kane unitary transformations 
T,j = el'^'T^ lie inside the even part Cl^{C) of C;„(C). 
Therefore, the Teo-Kane representation of Cm is just the 
spinor representation projectivized to Cm C SO{m). It 
follows that the projective image group of Teo-Kane repre- 
sentation of the even ribbon permutation group is Cm as an 
abstract group. Recall Clmi'C) is reducible into two irre- 
ducible sectors if m is even, and irreducible if m is odd. 
When TO is even, it is important to know how the relative 
action of the center Z{Gm) — Z2 of Gm on the two irre- 
ducible sectors. The center of Gm is generated by the ele- 
ment z = ((1 • • • 1), 1), whose corresponding Teo-Kane uni- 
tary transformation isU — 7172 • • • 7m up to an overall phase. 
As we show in Appendix H, the relative phase of z on the two 
sectors is always —1. 



VII. DISCUSSION 

We now review and discuss the main results derived in this 
paper. Using the topological classification of free fermion 
Hamiltonians^*"'^^, we considered a system of fermions in 3D 
which is allowed to have arbitrary superconducting order pa- 
rameter and arbitrary band structure; and is also allowed to 
develop any other possible symmetry-breaking order such as 
charge density-wave, etc. - i.e. we do not require that any 
symmetries are preserved. We argued in Section III that the 
space of possible gapped ground states of such a system is 
topologically equivalent to U{N)/0(N) for N large (for large 
N, the topology of U(A^)/0(A^) becomes independent of N). 
By extension, if we can spatially vary the superconducting or- 
der parameter and band structure at will with no regard to the 
energy cost, then there will be topologically stable point-like 
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defects classified by 7r2(t/(N)/0(N)) = Zj. 

This statement begs the question of whether one actually 
can vary the order parameter and band structure in order to 
create such defects. In a given system, the energy cost may 
simply be too high for the system to wind around U(A^)/0(A^) 
in going around such a defect. (This energy cost, which would 
include the condensation energy of various order parameters, 
is not taken into account in the free fermion problem.) If we 
create such defects, they may be so costly that it is energeti- 
cally favorable for them to simply unwind by closing the gap 
over large regions. (The energy cost associated with such an 
unwinding depends on the condensation energy of the order 
parameters involved, which is not included in the topological 
classification.) Thus, U(N)/0(N) is not the target space of an 
order parameter in the usual sense because the different points 
in U(N)/0(N) may not correspond to different ground states 
with the same energy. However, in Section II, we have given 
at least one concrete model of free fermions with no symme- 
tries in 3D in which the topological defects predicted by the 
general classification are present and stable. Furthermore, Teo 
and Kane^" have proposed several devices in which these de- 
fects are simply superconducting vortices at the boundary of 
a topological insulator 

In order to understand the quantum mechanics of these de- 
fects, it is important to first understand their quantum statis- 
tics. To do this, we analyzed the multi-defect configuration 
space; its fundamental group governs defect statistics. The 
configuration space of 2n point-like defects of a system with 
'order parameter' taking values in U(iV)/0(iV) is the space 
which we call i^2n- It can be understood as a fibration. The 
base space is X2n, the configuration space of 2n points (which 
we know has fundamental group S2n in dimension three and 
greater). Above each point in this base space there is a fiber 
Ai2n which is the space of maps from the ball B^ minu 
2n fixed points to 1]{N)/0(N) with winding number 1 about 
each of the 2n points. K2n is the total space of the fibra- 
tion. In Section V, we found that its fundamental group is 
7ri(i^2n) = Z X Z2 X S(Z^" X 5'2„). whcrc S(Z|" x Sm) is 
the subgroup of Z™ x S2n comprised of elements whose total 
parity in Zj" added to their parity in S2n is even. 

The fundamental group of the configuration space is the 
same as the group of equivalence classes of spacetime histo- 
ries of a system with 2n point-like defects. Since these differ- 
ent equivalence classes cannot be continuously deformed into 
each other, quantum mechanics allows us to assign them dif- 
ferent unitary matrices. These different unitaries form a repre- 
sentation of the fundamental group of the configuration space 
of the system. However, we found in Section VI that Teo and 
Kane's unitary transformations are not a linear representation 
of 7Ti{K2n), but n projective representation, which is to say 
that they represent 7Ti{K2n) only up to a phase. Equivalently, 
Teo and Kane's unitary transformations are an ordinary linear 
representation of a central extension of tti {K2n), as discussed 
in Section VI. 

This surprise lurks in a seemingly innocuous set of defect 
motions: those in which defects i and j are rotated by 2tt 
and the order parameter field surrounding them relaxes back 
to its initial configuration. This has the following effect on 



the Majorana fermion zero mode operators associated with the 
two defects: 



7i 



-It 



I3 



-Ij 



(56) 



One might initially expect that two such motions, one affect- 
ing defects i and j and the other affecting i and k, would 
commute since they simply multiply the operators involved 
by — signs. However, the unitary operator which generates 
(56) is^'': 



U^3^^tB^^^^ 



(57) 



Thus, the unitary operators U^^ and U do not commute; they 
anti -commute: 



jjij jjik _ _Jjik jjij ^ 



(58) 



However, as shown in Figure 12, the corresponding classical 
motions can be continuously deformed into each other. Thus, 
a linear representation of the fundamental group of the classi- 
cal configuration space would have these two operators com- 
muting. Instead, the quantum mechanics of this system in- 
volves a projective representation. 

Projective quasi-particle statistics were first proposed by 
Wilzcek-'", who suggested a projective representation of the 
permutation group in which generators (Tj and ak anti- 
commute for \j — fc| > 2, rather than commuting. Read^^ 
criticized this suggestion as being in conflict with locality. We 
can sharpen Read's criticism as follows. Suppose that one can 
perform the operation <ti by acting on a region of space, called 
A, containing particles 1 and 2, and one can perform 0-3 by 
acting on a region called B, containing particles 3 and 4, and 
suppose that regions A and B are disjoint. Consider the fol- 
lowing thought experiment: let Bob perform operation 0-3 at 
time t ~ and let Bob then repeat operation 0-3 at time t = 1. 
Let Alice prepare a spin in the state (1/V2)(| t) + l i)) at time 
t = — 1, and then let Alice perform the following sequence of 
operations. At time t = — e, for some small e, she performs 
the unitary operation | t) (^ | ® (Ti + 1 i) (i \<E)I, where / is the 
identity operation, leaving the particles alone. At time t = +e, 
she performs the unitary operation | t) (t | (8) / + | 4) (i | ® (Xi. 
Thus, if the spin is up, she performs ai at time t = — e, while 
if the spin is down, she does it at time t ~ +e. Finally, at time 
t = 2, Alice performs the operation cti again. One may then 
show that, due to the anti-commutation of ui and (T3, that the 
spin ends in the state (l/\/2)(| t) - | |)). However, if Bob 
had not performed any operations, the spin would have ended 
in the original state (l/\/2)(| t) + I i))- Thus, by perform- 
ing these interchange operations. Bob succeeds in transmitting 
information to a space-like separated region (if A and B are 
disjoint, and the time scale in the above thought experiment is 
sufficiently fast, then Alice and Bob are space-like separated 
throughout). 

Having seen this criticism, we can also see how Teo and 
Kane's construction evades it. The fundamental objects for 
Teo and Kane are not particles, but particles with ribbons at- 
tached. One may verify that, in every case where operations 
in Teo and Kane's construction anti-commute, the two opera- 
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tions do not act on spatially disjoint regions due to the attached 
ribbons. That is, the interchange of particles also requires a 
rearrangement of the order parameter field. 

While this argument explains why a projective representa- 
tion does not violate causality, it does not really explain why 
a projective representation actually occurs in this system. Per- 
haps one clue is the fact that the hedgehogs have long-ranged 
interactions in any concrete model. Even in the 'best-case sce- 
nario', in which the underlying Hamiltonian of the system is 
U(A^)/0(iV)-invariant, there will be a linearly-diverging gra- 
dient energy for an isolated hedgehog configuration. Thus, 
there will be a linear long-ranged force between hedgehogs. 
Consequently, one might adopt the point of view that, as a re- 
sult of these long-ranged interactions, the overall phase asso- 
ciated with a motion of the hedgehogs is not a purely topo- 
logical quantity (but will, instead depend on details of the 
motion) and, therefore, need not faithfully represent the un- 
derlying fundamental group. As one motion is continuously 
deformed into another in Figure 12, the phase of the wave- 
function varies continuously from +1 to — 1 as the order pa- 
rameter evolution is deformed. It is helpful to compare this 
to another example of a projective representation: a charged 
particle in a magnetic field B. Although the system is in- 
variant under the Abelian group of translations, the quantum 
mechanics of the system is governed by a non-Abelian pro- 
jective representation of this group (which may be viewed as 
a linear representation of the 'magnetic translation group'). 
A translation by a in the x-direction, followed by a transla- 
tion by b in the y-direction differs in its action on the wave- 
function from the same translations in the opposite order by a 
phase abB/^o equal to the magnetic flux through the area ab 
in units of the flux quantum $0. If we continuously deform 
these two sequences of translations into each other, the phase 
of the wavefunction varies continuously. For any trajectory 
along this one-parameter family of trajectories (or 'movie of 
movies'), the resulting phase of the wavefunction is given by 
the magnetic flux enclosed by the composition of this trajec- 
tory and the inverse of the initial one. In our model of non- 
Abelian projective statistics, the phase changes continuously 
in the same way, but as a result of the evolution of the order 
parameter away from the defects, rather than as a result of a 
magnetic field. 

An obvious question presents itself: is there a related theory 
in which the hedgehogs are no longer confined? Equivalently, 
could 3D objects with non-Abelian ribbon permutation statis- 
tics ever be the weakly-coupled low-energy quasiparticles of 
a system? The most straightforward route will not work: if we 
had a U(A^)-invariant system and tried to gauge it to eliminate 
the linear confining force between hedgehogs, we would find 
that the theory is sick due to the chiral anomaly. If we dou- 
bled the number of fermions in order to eliminate the anomaly, 
there would be two Majorana modes in the core of each hedge- 
hog, and their energies could be split away from zero by a 
local interaction. This is not surprising since ribbon permu- 
tation statistics would violate locality if the hedgehogs were 
truly decoupled (or had exponentially-decaying interactions). 
On the other hand, if the hedgehogs were to interact through 
a Coulomb interaction (or, perhaps, some other power-law). 



they would be neither decoupled nor confined, thereby satis- 
fying the requirements that they satisfy locality and are low- 
energy particle-like excitations of the system. Elsewhere, we 
will describe a model which realizes this scenario*"" 

The non-Abelian projective statistics studied in the 3D class 
with no symmetry can be generalized to arbitrary dimension. 
As shown in Eq. (41), the classification of topological defects 
is independent of the spatial dimension. Thus, in any dimen- 
sion d with no symmetry (p ~ 0), point-like (D ~ 0) topolog- 
ical defects are classified by TTa{Rp-D+i) = t^o{Ri) = ^2- 
Moreover, it can be proved that analogous topological defects 
in different dimensions not only carry the same topological 
quantum number, but also have the same statistics. In Sec. 
V, we have defined the configuration space A^2n which is 
the space of maps from i?'^ \ lJi"i ^f ^o R7 = U/O, with 
specific boundary conditions. Now if we consider point de- 
fect in the class with no symmetry in 4D, the configuration 
space A^2n'' is defined by maps from _B^ \ IJ^"]^ Bf to the 
classifying space i?g = Sp/U. Noticing that i?^ \ IJ^"]^ Bf 
is homotopy equivalent to the suspension of B'^ \ lJi=i Bf, 
we obtain that M'^n'^ i^ equivalent to the space of maps 
from B^ \ U,^"i Bf -^ fl {Sp/U), where fl (Sp/U) is the 
loop space of Sp/U. Since il [Sp/U) ~ U/O, we obtain 
A^2n* — A^2n- Thus wc havc proved that the configuration 
space M.2n is independent of the spatial dimension d. On the 
other hand, the fundamental group of the configuration space 
X2ri of 2n distinct points in B'^ is independent of d as long 
as d > 2. Consequently, the space K2n defined by the fibra- 
tion A^2n ^ K2n — > X2n is also topologically independent 
of spatial dimension d for d > 2. Thus, the proof we did 
for d — i applies to generic dimension, and non-Abelian pro- 
jective statistics exist in any spatial dimension d > 3 for point 
defects in the no-symmetry class. A similar analysis applies to 
extended defects with dimension I? > 0. When the spatial di- 
mension is increased by 1 and the symmetry class remains the 
same, the classifying space is always changed from i?2+p-d 
to R2+p-d-i — ^^^{R2+p-d)- Consequently, at least for 
simple defects with the topology of S^ , the statistics is inde- 
pendent of spatial dimension d as long as d is large enough. 
For point defects, the "lower critical dimension" is d = 3, 
while for line defects, i.e. D = 1, the "lower critical dimen- 
sion" is at least d — A since in rf = 3 we can have braiding 
between loops. 



Appendix A: Topological Classification of Hamlltonlans with 
Charge Conservation but without Time-Reversal Symmetry 



For the sake of completeness, in this appendix we give the 
topological classification of free fermion Hamiltonians with 
charge conservation but not time-reversal symmetry. We be- 
gin with the zero-dimensional case. The charge conservation 
condition is most transparent when the Hamiltonian is writ- 
ten in terms of complex fermion operators, c^, c|, rather than 
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Majorana fermions: 



H 



j:h, 



ic]c. 



(Al) 



«j 



The indices z, j range from 1 to N, the number of bands. So 
long as there is one c and one c^ in each term, the Hamihonian 
will conserve charge. We make no further assumption about 
the Hamihonian. As in our discussion in Section III, we are 
only interested in the topology of the space of such Hamil- 
tonians, so two Hamiltonians are considered to be equivalent 
if they can be continuously deformed into each other without 
closing the gap. Thus, we can flatten the spectrum and assume 
that H and, therefore, hij only has eigenvalues ±1; any other 
gapped Hamihonian can be continuously deformed so that it 
satisfies this condition. Then, hij can be written in the form; 



/l 



h^U^ 



\ 



-1 



V 



u 



(A2) 



-1/ 



with k diagonal entries equal to +1 and N—k diagonal entries 
equal to —1 for some k. The space of such matrices hij is 
equal to the space of matrices U , namely U( A^), modulo those 
matrices U which commute with the diagonal matrix above, 
namely U(fc)xU(A^ — k). Thus, the space of matrices hij is 
topologically equivalent to 



N 



M2N = U U(iV)/(U(fc) X \]{N - k)) 



(A3) 



k=0 



For N sufficiently large, the topology of these spaces cannot 
depend on k, so we write this simply as 'Lx\J{N)/{\J{N /2)x 
U(7V/2)). 

Now suppose that we have a ID system. As in section III, 
we will consider the Dirac equation which approximates the 
Hamihonian in the vicinity of the points in the Brillouin zone 
at which it becomes small. We write: 



iJ = V^(i7i5i+M)V' 



(A4) 



Here, we have suppressed the band index on the fermion op- 
erators ijj^ and if). Both 7 and M are N x N matrices; 71 
satisfies tr(7i) = and 7J = 1. Then, the different fermion 
operators will have the same gap if 



{7i,M} = 0. 



M^ 



(A5) 



Viewed as a linear operator on C , 71 defines two N/2- 
dimensional subspaces, its eigenvalue +1 eigenspace and its 
eigenvalue —1 eigenspace. Let us call them X+ and X_. 
Since {71, M} == 0, M is a map from X+ to X_. Thus, M is 
an isometry between two copies of C^^^ and can be viewed 
as an element of U(iV/2). This can be made more concrete by 



taking, without loss of generality. 



7i = 



I 

-I 



(A6) 



where I is the N/2 x N/2 identity matrix. Then, M can be 
any matrix of the form 



M 



W 
U 



(A7) 



fovU e\J{N/2). 

Now, consider a 2D system. The Dirac equation takes the 
form: 



H ^xP^{ijjdj +M)tP 



(A8) 



with j = 1,2. As before, 71 defines two A^/2-dimensional 
subspaces, X+ and X_. Then i"j2M commutes with 71 and 
squares to 1. Thus, it divides X^ into two spaces, X\_ and 
X"^ which are the +1 and —1 eigenspaces for J72M (and it 
similarly divides for XJ). Thus, a choice of M is a choice of 
linear subspace X\_ of X+ or, in other words, 

N 

M € IJ U(iV/2)/(U(fc) X U(iV/2 - k)) 

k=0 

= Z X U(A^/2)/(U(A^/4) X U(A^/4)) (A9) 

Continuing in this way, we derive Table III analogous to Ta- 
ble I but for charge-conserving systems without time-reversal 
symmetry. Table III has period-2 as the dimension is increased 
whereas Table I had period- 8. 



The two rows in Table III correspond to the presence or ab- 
sence of an additional unitary symmetry which requires the 
eigenvalues of the Hamihonian to come in pairs with equal 
and opposite energy ±Ei, i.e. a sublattice symmetry, as dis- 
cussed in Sec. Ill F. Consider a zero-dimensional system de- 
scribed by the Hamihonian (Al) with such a symmetry. Then 
there is a unitary matrix U satisfying 



{[/, h}^0 



U^ = 1 



(AlO) 



We can now apply the same logic to U and M that we ap- 
plied to 7i and M when considering one-dimensional sys- 
tems. Thus, M e U{N/2). More generally, it is clear that 
U plays the role of an extra 7-matrix, so the row of Table III 
for systems with sublattice symmetry is just shifted by one 
from the row of the table for systems without it. 

In order to compute the homotopy groups of these classify- 
ing spaces, it is useful to follow the logic which we employed 
in Sections IIID,F. We approximate the loop space of U{N) 



by minimal geodesies from I to —I: L{X) 



„i\A 



where A 



is a Hermitian matrix satisfying A = 1 so that L{Q) =1 and 
L{n) = —I. Any such matrix A divides C^ into +1 and —1 
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dim.: 

no symm. 

sublattice symm. 



Z X 





U(Af/2)xU(Af/2) 
\5{N/2) 



Z X 



1 

U(iV/2) 

U(N/2) 
U(JV/4)xU(Ar/4) 



Z X 



U(JV/2) 
U(Ar/4)xU(Af/4) 

U(iV/4) 



Z X 



\]{N/2) 

U(JV/4) 
U(JV/8)xU(Ar/8) 



TABLE III: The period-2 (in both dimension and number of symmetries) table of classifying spaces for charge-conserving free fermion 
Hamiltonians without time-reversal symmetry. The systems have A'^ complex fermon fields in dimensions d = 0,l,2,3,... with no additional 
symmetries (apart from charge-conservation) or with an additional symmetry, such as a sublattice symmetry. Moving p steps to the right and 
p steps down leads to the same classifying space (but for 1/2*" as many fermion fields), which is a reflection of Bott periodicity, as explained 
in the text. The number of disconnected components of any such classifying space - i.e. the number of different phases in that symmetry class 
and dimension ~ is given by the corresponding ttq, which may be found in Eq. A13. Higher homotopy groups, which classify defects, can be 
computed using Eq. A12. 



Symmetry classes 


Physical realizations 


d=l 


d = 2 


d = 3 


A 


Generic ins. 





Quantum Hall (GaAs, etc.) 





AIII 


Bipartite ins. 


Carbon nanotube 





Z 



TABLE IV: Topological periodic table in physical dimensions 1,2,3 of the complex classes. A and AIII in the first column labels the 
symmetry classes in the Zirnbauer classification.^"*''^. The second column are the requirements to the physical systems which can realize 
the corresponding symmetry classes. The three columns d = 1, 2, 3 list the topological states in the spatial dimensions 1, 2, 3 respectively. 
means the topological classification is trivial. The red label Z stands for the topological state which have not been realized in realistic materials. 



eigenspaces or, in other words, it can be written in the form 



A = C/^ 



-I 



u 



(All) 



N-k 



where I,„ is the m x m identity matrix. Therefore A E 

^ ^ ufjv/2')xufJV/2') - Thus, the loop space of \J{N) can be 
approximated by Z x yji^N/2)xv{N/2) ■ ^hus, 7rfc(U(iV)) = 
^fe_i(Z X l]{N)/{\J{N/2) X V{N/2))). 

To compute homotopy groups of Z x U(A^)/(U(A^/2) x 
U(A^/2)), consider the homotopy long exact sequence associ- 
ated to the fibration 

U(iV/2) -^ U{N) -^ V{N)/U{N/2). 

It tells us that 

. . .^,_i(U(7V/2)) ^ ^,(U(iV)) ^ 7r,(U(A^)/U(7V/2)) .... 

If we are in the stable limit, i.e. for N sufficiently large, then 

^,_i(U(7V/2)) = 7r,(U(A^)), so ^,(U(iV)/U(7V/2)) = 0. 
Now consider the fibration 

U(iV/2) -^ \]{N)/V{N/2) -^ U{N)/{\]{N/2) x U(iV/2)). 

Then, since 7ri{\J{N)/\J{N/2)) = 0, the associated homo- 
topy exact sequence is simply 

^ n,{\]{N)/{U{N/2)x\]{N/2))) ^ 7r,_i(U(iV/2)) ^ 

so7r,(U(7V)/(U(7V/2) X U(iV/2))) = 7r,_i(U(iV/2)). 

Combining these two relations, we have Bott periodicity for 
the unitary group: 

TTk{U{N)) = TTk-iiZ X U{N)/{U{N/2) X U(7V/2))) 

= TTk-2{V{N/2)) (A12) 



By inspection, we see that 

7ro(U(7V)) = 
7ro(Z X U(iV)/(U(iV/2) x U(iV/2))) = Z (A13) 

Combining this with Eq. A 12, we learn that the even homo- 
topy groups of V{N) vanish and the odd ones are Z; the re- 
verse holds for U(A^)/(U(A^/2) X U(A^/2))). 

These homotopy groups can be used to classify the dif- 
ferent possible phases and topological defects of systems in 
the symmetry classes/dimensions in Table III. Similar to Ta- 
ble II in Sec. IIIF, we list in Table IV the nontrivial topo- 
logical physical systems in the two symmetry classes. In 
Ref. 26, these were called A and AIII, respectively, following 
the corresponding classification in random matrix theory^"^'^^ . 
In A class (with charge conservation and without sublattice 
symmetry or time-reversal symmetry) a topological nontriv- 
ial state is classified by integer in 2d, which is the famous 
quantum Hall system. In the AIII class with sublattice sym- 
metry, in Id the carbon nanotube or graphene ribbon can be 
considered as an example with nontrivial edge states. Com- 
pared to the requirement of this symmetry class, carbon nan- 
otube has too high symmetry because time-reversal symmetry 
is present. An orbital magnetic field can be coupled minimally 
to the carbon nanotube to get an example system with exactly 
the required symmetry of AIII class. In 3d there should be an 
integer classification of the AIII class but no realistic topolog- 
ical nontrivial material is known yet. 



Appendix B: Understanding the Role of Spatial Dimension in 
the Classification of Free Fermion Systems 

In Sections III and A, we have seen that increasing the 
spatial dimension moves a system through the progression 
of classifying spaces oppositely to adding symmetries which 
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square to — 1 (but in the same direction as adding symmetries 
which square to 1). The classifying space for d-dimensional 
systems is the loop space of the classifying space for d + 1- 
dimensional systems with the same symmetries. Symboli- 
cally, 



Thus, for lattice-translational symmetry, we may employ 



R2 



0—d 



n{R: 



■2+p-(d+l) 



(Bl) 



where fl denotes the loop space and, as in Section III, p is the 
number of symmetries and d is the spatial dimension. Thus, 
in using Tables I and III, moving one step to the right, which 
increases the dimension and one step down, which increases 
the number of symmetries, leaves the classifying space the 
same. 

We explained this in Sections III and A by expanding the 
Hamiltonian about the point(s) in the Brillouin zone where the 
gap is minimum where it has the form of the Dirac Hamilto- 
nian. Analyzing the space of such Hamiltonians then amounts 
to analyzing the possible mass terms, namely the space of 
matrices which square to —1 and anticommute with the 7- 
matrices and anti-unitary symmetry generators. One might 
worry that this analysis is not completely general since it de- 
pends on our ability to expand the Hamiltonian in the form of 
the Dirac Hamiltonian, and this seems non-generic. However, 
Kitaev's texture theorem (unpublished) states that a general 
gapped free fermion H canonically deforms to a Dirac Hamil- 
tonian as above without closing the gap, so only these need to 
be considered. 

Another more algebraic approach involves the suspen- 
sion isomorphism in KR theory. In complex K-theory: 

KV+'^{X A S"^) ~ KU^-iX), however, for KR[A] with 
'_{xi,...,Xd) == (-.Ti,...,-a:;d), one gets: KR^-'^{X A 
S"^) ~ KO"{X). (There is a similar isomophism in twisted 
K-theory where (^)^ = —1.) In (translation invariant) parti- 
cle/hole symmetric systems 



-^^jUA^,-x,7x, = -^^lxAx,~-x,ll^ 



= ^^ /_^lxj^X,-Xjlx, = — ^ /_^lxi^X,~XjlXj7 

so reversing a space coordinate introduces a minus sign, it can 
be argued via a high frequency cutoff that momentum space 

directions correspond not to M. but actually to S . 

If translational symmetry is broken to a lattice Z'^ sym- 
metry, then momentum space becomes a d— torus T"^ and it 
should replace S"^ in the l.h.s. of the two isomorphisms above. 
Fortunately, all K-theories are generalized cohomology theo- 
ries and thus depend only on stable homotopy type. Tori have 
extremely simple stable types: 



E^' 



(B2) 
where ^ denotes suspension. The same formula holds with 
(^) above all spaces. 



KBJ'-'^{X AT'^) 



KO (X). 



""''^^X) ffi ifO"-''(X). (B3) 



®/d-j KO 



In this appendix, we want to sketch here a third way of un- 
derstanding "spatial dimension = de-loop" that is more geo- 
metric than either of the alternatives above. 

We warn the reader that this section is only a "sketch", so 
there may be tricky analytic details regarding the precise def- 
inition of the controlled spaces which we have over-looked. 
The math literature^^ in controlled K-theory is constructed in 
a less naive context. We are not sure if there is something 
essential we are missing. 

Let Oz be a limitingly large orthogonal group with integer 
(Z) control and a topology (like compact/open) where sliding 
a disturbance off to ±00 converges to the identity. That is, Oz 
is the set of infinite matrices O/^j with a finite-to-one map 
(decapitahzing) J ^ j <E Z such that \Oij\ < 6-™"=*'!*-^! 
and O'j jOi^K = 6.J.K with the following identifications: 
Any matrix O/.j with O/q^j = ^lo.J is equivalent to the ma- 
trix with one less row and column obtained by removing the 
/* row and /* column. 

By changing I, the number of distinct J map that map to a 
given i, we can choose any desired constant in the exponential 
decay cxp(— constji — j|). We could also have chosen to study 
the class of infinite matrices with O/.j = for |i — j| > 1; 
in this case, the matrix O/j is a banded matrix, which is non- 
zero only within distance / of the main diagonal. We refer to 
the case of exponential decay as "soft control", while we refer 
to the case that O/^j = for |j — j| > 1 as "strict control". 

The key claim in this "third" approach is 



Claim. Oz ~ '^/q ^qxZ = BOxZ(z 
homotopy type.) 



denotes weak 



We take for BO x Z the well known model of pairs, con- 
sidered as formal differences, of transverse (not necessarily 
spanning) subspaces (L_|_ , L_ ) of a d-dimensional real vector 
space R'^ (with d finite and a limit taken as d -^ 00.) The 
Z-coordinate above is the index dim_L+ — dimL_. 

Here are the maps: choose any cut in Z, say 11+ projects to 
J with j > and n_ projects to J with j < 0. Set 

/ : Oz ^ BO X Z, 

/(O) = (kcr(n+ o OI5J, ker(n_ o 0\s_)), 

where 5'+(_) is the space of basis elements mapping to 
[0, L] ([— £, —1]) where I ^ the decay constant for Oz, and 
both kernels are regarded as subsets of the span M."^ of ele- 
ments mapping to [—L, L]. Given (i+, L_) G BO x Z, de- 
fine g{L+,L^) e Oz as follows. Think of L± C W^ with 
P denoting the perpendicular subspace to span(i+, L_). The 
infinite vector space on which g{L+, L_) acts is spanned by 
a copy of W^ for each i € Z. O — g{L^,L_) acts as the 
identity on each P x i, translation by — 1 e Z on each _L+ x i 
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and translation by +1 G Z on each L_ x i. Note: the use of 
the kernel above is correct in the case of strict control. In the 
case of soft control, we should replace the kernel by the pro- 
jector onto right singular vectors of 11+ o 0\s+ with singular 
value less than or equal to some quantity which is exponen- 
tially small in L. 

It is immediate that fog — Idsoxz- On the other hand, go 
f seems to turn a general O G Oz into a very special form (no 
rotation at all, just various subspaces sliding left and right.) 
However, we claim that there is a deformation retraction of Oz 
to maps which take this simple form: r : Oz — >■ {slides} C 
Oz- First, think about a fiberwise O (no sliding) which is X 
over some i E Z and Id over other j ^ i E Z. Note that 



X 







Id 





-X 




Id 



via the canonical homotopy 



TTi 



X 
-X 



TTt 



Id 
Id 



When d is even (which we may assume) there are further 
canonical homotopies so that 



Id 
Id 



<Jz 
az 












(TZ 


^ 


Id 


(Tz 







Now a "Eilenberg swindle" runs X (canonically) towards 

+00. 



Id X Id 



lA X -X X -X 



Id Id Id Id ■■■ 

^ I (i + !)(« + 2)(« + 3)- •• 



4). 



Reusing this trick, any fiber-wise, though now it is better to 
slide from the center out to both ±00, X can be canonically 
connected to Id G Oz- Similarly, one can (canonically) de- 
form the general O G Oz to a sum of three pieces — fixed, 
left-sliding, right-sliding — and this is the deformation retrac- 
tion r written above. 

We now promote the (weak) homotopy equivalence Oz — 
BO X Z to an entire table where the rows are all (weak) ho- 
motopy equivalences. 



^^VUz^ O1 



O7, 



—2 



Oz 



Vo, 



z+z 



xOz 
Oz 



a codimension 1 hyperplane, in this case the line y = —\, 
and finds a pair of locally finite Z-controUed kernels which 



define an element in 



Oz 



yu,. 



% 

ii BO xZ 

:^3 O 



O 



Vu 



Moving diagonally up and left, say from ~i to ~2 seems 
to be easy: one cuts the control space not at a point but along 



O 



/ Cz X Cz 



However, there is 



one subtle point: the matrices 11+ and O are both controlled; 
however, this does not imply that the projector onto the kernel 
of n_|_ o 0\sj^ is also controlled. If we knew that the singular 
values of 11+ o 0\sj^ had a spectral gap separating the zero 
singular values (or, in the case of soft control, the singular 
values which are exponentially small in L) from the rest of 
the spectrum, then we could show that the projector onto the 
kernel had soft control, with decay constant set by the gap, but 
in absence of knowledge of a gap, this seems to be a difficult 
technical step. 

Moving down the array, say from ~i to ~3, requires "loop- 



O, 



i.e. that the 



ing" ~i. We must check that Vl{Oz) — '^'Vjj , 
first step in Bott's ladder of eight rungs holds with Z-controls. 
To do this, one must do Lie theory in Oz- find the shortest 
geodesic arcs from Id to —Id in Oz, understand which direc- 
tions they pick out in the Lie algebra Oz and that their mid- 
points are Z-controlled complex structures, and finally make 
the analogs of Bott's index calculations. The initial step is to 
consider a representative for a "Z-controlled complex struc- 
ture" (where coordinates are consecutive in Z): 



J = 0^ 



1 




-1 






1 




-1 



o, 



OgOz 



and see that the 1 -parameter subgroup {e* '"^'^j C Oz retains 
Z-control. This may be checked by writing 



log J 



-J, 



since J^ ~ — /. This yields Z-controls on log J and hence on 
the entire 1 -parameter subgroup. 

One interesting feature of this approach is that we are able 
to classify different symmetry classes of controlled unitaries 
on a line. The controlled unitaries Uz are classified by an inte- 
ger. This integer is precisely the "flow" described by Kitaev^^. 
If the controlled unitary is assumed to be symmetric, it is in 
the orthogonal group, and again Oz is classified by an integer. 
However, if the unitary is chosen to be self-dual, in this case 

^^Sv ^^^ ^ ^2 classification. Ryu et. al.^*" note (see Table 
I or Table IV of that paper) that each of the 10 different sym- 
metric spaces can be obtained by considering the exponential 
of Hamiltonians lying in the 10 different symmetry classes; 
for example, given a Hermitian matrix B, the exponential 
cxp{iHt) is a unitary matrix, while if B is anti-symmetric and 
Hermitian, the exponential cxp{iHt) is an orthogonal matrix. 
However, the classification of unitaries here implies that there 
can be an obstruction to writing a controlled unitary as an ex- 
ponential exp{iHt) with H of the correct symmetry class and 
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H also controlled. 



where a — {ax, o-y,(Jz) are Pauli matrices and the spinor 



Appendix C: Phase Symmetry Unbroken 



(D2) 



If the phase symmetry of the order parameter is unbroken 
(i.e. if we allow gradients in the overall phase), then the ef- 
fective sigma model target is U/O, not U/O. This means 
that at every stage of the calculation, an additional circle ^rjl 
or K{'L,1) product factor must be added to the target space. 
This leads to no change for 7r2(y) or 7ri(_R) (essentially be- 
cause 7r2(QTj4) = 0), however iriiY) and tto{R) pick up an 
additional integer which is the winding number around this 
circle. Thus, (CI) below 






l\\ 



m 



^i" >^ S2n ^2 



(CI) 



is replaced by (C2): 



a 



^^(ynew) ^^^(^new) V7ri(if2^") 




7ri(y"^") ^7ro(i?' 

anew 

(Z2" X Z2") X S2n Z^" X Z2 

y2n 



(C2) 



where S2n acts on Z by permuting its standard generating 
set. 

The new Z^" factor in 7ro(-R"''™) represents rolling the or- 
der parameter around the phase circle as one moves radially 
into the hedgehog. Similarly, the Z^" in tti{Y"™) comes 
from wrapping the whole hedgehog around the phase circle 
as the parameter is advanced. Clearly, the two Z^"'s cancel: 
ker(a5'™) == ker(9i). Thus, 

TTliK2n) = 7ri(^2«) X Z, 

where the extra factor Z comes from a base point in B'^ wraps 
around the additional circle jjij4 . 



Appendix D: Hopf Map: an example of boundary maps in the 
homotopy exact sequnce 

The Hopf fibration is an example of a fibration (of a fiber 
bundle, in fact) which arises in several contexts in physics. It 
consists of a map (the Hopf map) from S^ to S*^ such that 
the pre-image of each point in S*^ is an S^. This is familiar 
from the representation of a vector (e.g. the Bloch sphere) by 
a spinor (e.g. a two-state system): 



satisfies |zip + |z2|^ = 1 so that n is a unit vector. Due to 
this normalization condition, the spinor z lives in S^ while 
the vector n lives in 5*^. Transforming z ^ e*^ z leaves n 
invariant. Thus, the map (Dl) maps circles in S^ to points in 
S^, as depicted below: 

t/(l) phase = s-i ^ 5-3 = normahzed states in C^ 
5*2 = Bloch sphere 

One context, quite closely-related to the present discussion, 
in which the Hopf map arises is the 0(3) non-linear <t model 
in 2 + 1-dimensions. The order parameter field n(a;) is a map 
from S^ to S^ (assuming that n(a:) is constant at infinity, so 
that spacetime can be compactified into S^). The topological 
term which can be added to such a model^^ is: 



S 



H 



1 

'2^ 



d^xdref'-'^A^F.x 



where Fi^x = d„Ax — d\Ai, and 



e^^'^'-d^Ax 



1 

8^ 



e^"^^!! • d^n X dxn 



(D3) 



(D4) 



This term gives the linking number of skyrmion trajectories 
in the 0(3) non-linear a model. To see this, suppose that 
n(oo) — (0,0,1). The pre-image of (0, 0, —1) (which we 
may view as the centers of skyrmions) is a set of closed cir- 
cles, and the Hopf term gives their linking number. 

We can use the Hopf map to explain boundary maps in the 
homotopy exact sequence: tt2{S'^) — > 7ri(5'^): 



211 



211 



211 




Z' (7Z 



(Dl) 



which, we see, must be an isomorphism to be consistent with 
TTi {S^) ^ 0. But how is this boundary map defined? We think 
of a 2-sphere S"^ as a loop of circles S^, as depicted in Fig. 14. 
Using the (defining) property of a fibration (see any book on 
algebraic topology), we may lift this loop of circles S^ from 
the base S^ to the total space S^ . The lifting is perfectly con- 
tinuous, but it does not give us a sphere. Instead, it gives us 
a surface with boundary. This boundary is an S^. What hap- 
pens is that, at the very end of the loop of circles 5*^, the lifting 
"explodes" the final circle 5*^, whose map to the Bloch sphere 
is shrinking to a point, into a circle in 5*^, as pictured in figure 
15. In other words, an element of tt2{S'^) is lifted to a surface 
in S^ whose boundary is a circle. This circle can wind around 
the fiber 5*^, thus defining an element of 7ri(S'^). This map 
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base 5* = Bloch sphere 



FIG. 14: The circle S^ starts small, gets large, and becomes small 
again at the end of the "loop," that is, after it has swung all the way 
around the sphere S^. 



from an element of 7r2(5'^) to an element of 7ri(S'^) is called 
a boundary map since the element of 7ri(S'^) is the boundary 
of the hft to S^ of the corresponding element of 7r2(S'^). 



another fiber 



r 



S"^ -< 



V. 




fiber over base pt. — 
lift at end of family 



base pt. = 
lift at beginning 
of family 



S'^ \ A^ is already contractible, to make S^ IJ^^ 2 ^^ ^'^^' 
tractible, we must remove a diameter S C D^ whose end- 
points are glued to N. Making our maps / transverse to 
S C D^ (and not just N C S"^ as would be ordinarily done 
in the PTC) we see that the ribbons (and their framed cobor- 
disms) are now subject to a "singularity" at the /-preimages 
of the origin O C D^. The singularity is very mild: it is sim- 
ply a point on the ribbon where the orientation (induced by 
comparing orientations on M^ and 5^) longitudinally along 
the ribbon ro changes direction. The upshot is that the gener- 
ators for 7ro(-R2n) now consist of ribbons with no longitudinal 
arrows assigning a direction along the ribbon. (These gener- 
ating ribbons are, of course, subject to the boundary condition 
that there are exactly 2n arc endpoints meeting each hedge- 
hog exactly once at its marked point f^^(N).) Closed loops 
of ribbons may also arise. Because all maps from a 3-manifold 

deform into these lowest cells, 5^ IJ^^ ^ ^^ C U/O, all gen- 
erators of 7ro(i?2n) are of this form. 

Finding the relations defining 7ro(i?2n) requires looking at 
maps of a four manifold {Q x /, d{Q x /)) into U/O. To 
study these by a variant of the PTC, we would need to take 

into account all cells up to dimension 4, sk^{U /O), the "4- 
skeleton." Similarly, a full analysis of tti (R2n) would require 

looking at sfc'^(C//0). 



FactE.l. 7r3(52U,,,2 
Hopfmap h:S^^ S^. 



D^ 



and is generated by the 



FIG. 15: The lifts of circles S^ are depicted 'upstairs' with respect to 
the usual picture of Hopf circles filling K^. The lift at the end of the 
family is the biggest ellipse in this picture, not a single point, unlike 
the lift at the beginning of the family. Thus, the loop of circles no 
longer closes when lifted. 



Proof. The Whitehead exact sequence reads: 



Hi{X) ^r(7r2(X)) ^^aW "^'^^^^^ > ffgW 



Appendix E: Computation of tto (Rin), t^i {R2n) using the 
"Postnikov tower" for U/O 

In this appendix, we compute the homotopy groups 

7I'o(^2n),7ri(i?2n)- 

Recall that the homotopy classes of maps from a 3-manifold 
M^ to S^ can be understood in terms of the ribbons which 
are the (framed) inverse image of the north pole N C S'^. 
Such ribbons serve as generators for the path components of 
Maps (A'/'' -^ S*^) and the relations are "framed cobordisms," 
that is, imbedded surfaces (5;ro,ri) C (M^ x [0,1]; M x 
0, M X 1) where S* is a surface "cobording" between initial 
and final ribbons ro and ri. S, like the ribbons themselves, 
has a normal framing which restricts to the oriented ribbon 
directions of both ro and ri . This is a special case of what, in 
topology, is called the Pontryagin-Thom construction (PTC) 
(see Appendix F for more details). 

Of course, there is more to U/O than the bottom S"^. Let's 
look next at the influence of the 3-cell D^ on Ron. Whereas 



where F is Whitehead's quadratic functor. It is known that 

F(Z2) = Z4. If X = 52 U^^g 2 -D^ the outer integral homol- 
ogy groups vanish, so ir^iX) = Z4 as claimed. 

To see that the generator g E Tr3{X) deforms to S'^, make 
transverse to the origin O C D^. Since D^ is attached 
to S'^ with positive degree, the signed sum of inverse im- 
ages g^^(O) = 0. These may be paired and then removed 
by building a framed cobordism (see PTC appendix) from 
g^^(O) to inside S^ x [0, 1]. The map g' associated to the 
5"^ X 1 level avoids O C D^ and so may be radially deformed 
into 5^. 

n 



Consequence E.2. .sk%U/0) = S^ Urf.,2 ^' U2HopfD^ 



Proof. Consider the commutative diagram: 
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^3 ^ SU{2) ■ 



SO 
fiber 

Y 

-^su 



Hopf map 



coset 



coset 



S 



2 ^ 



SU{2) 



SO(2) 

Applying the tt^ functor, we obtain: 



su 
so 




epimorphism 



From the homotopy exact sequence of 
SO^SU 



i; =u/o 



we can see that 

7r3(50)^^7r3(S't/)- 
l\\ l\\ 



-7i"3( 



SO' 



-^772(80) 

1 



and we conclude that 5 is an epimorphism. Thus, the Hopf 
map, which is order 4 in X, needs a 4-cell to be attached in 

order to make it order 2 in U/O. O 

We will not need this, but for the curious, sk^{U/0) re- 
quires one additional 5-cell which maps degree 2 over the 4- 
cell. 

As remarked, studying maps into cell structures becomes 
laborious. Homotopy theory works best when studying maps 
out o/cell structures and into fibrations. So, let us introduce 

the 2-stage Postnikov tower for U/O. 

The notation K{'!t, n) is used for any connected space (gen- 
erally infinite dimensional) that has only a single nontrivial 
homotopy group, 7r„(iir(7r, n)) = it. It is easy to show that 
the homotopy type of K{'k, n) is unique. At the next stage 
of complexity come spaces Z with two nontrivial homotopy 
groups, say tt„i{Z) = A and 7Tn{Z) = B, m < n. Such 
spaces (for simplicity we now assume m > 2) have a homo- 
topy type determined by the total space T of a fibration called 
a 2-stage Postnikov Tower. 



K{n,B) 
(fiber) 



-> T 



K{m,A) (base) 

To completely specify how the fibration twists and the 
homotopy type of T we need to specify the "fc" invariant 
fee H'^+^{K{m,A);B). 

The space K{m, A) classifies cohomology in the sense that 
for any space S, the homotopy classes of maps [S, K{m, A)] 
naturally biject with H"^{S;A). (There is a "fundamental" 
class L e H"^{K{m, A); A) so that for / : 5 ^ K{m, A) we 
associate /*t G H"^{S; A).) Thus, the k invariant is really a 
map: 

k : K{m,A) -^ K{n + 1,B) 

and T is the pullback of the path loop fibration over K{n + 

l,B): 



nK{n + 1, B) = K{n, B) > pt. 



K{n,B) 



K{n+l,B) 



K{n,B) 



K{m,A)- 



k 



-^pt. 



■Kin + 1,B) 



where we have denoted the contractible space of paths in 
K{n + 1,B) starting at its base point by its homotopy model, 
a single point, pt. 

For spaces with many homotopy groups, the Postnikov 
tower can be continued iteratively. 

In our case, TTi{U /O) — 0,Z2,Z2,0 for i == 1, 2, 3, 4, so a 
2-stage tower is adequate for computing TTi{R) for i = 0, 1. 

Thus, our model T for U/O is now a fibration: 



if(Z2,3)- 



-> T 



if(Z2,2) 
Similar to the previously computed cell structure of U/O, 
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the cells of K {'1^2,2) through dimension < 3 are the same. 
But for K{Z2,2), the 4-cell kills the Hopf map and the 5-cell 
is now degree 4 over the 4-cell. The cell complex for K{Z2, 2) 
begins as follows: 



to R but with maps to T replaced with maps to its base space 
K{Z2,2). Similarly, RS is Maps(Q -^ K{Z2,3)) with dQ 
mapping to the base point of K{Z2, 3). We have from the 
homotopy exact sequence: 



Cs- 


>C^- 




-^C3- 


>C2 


l\\ 
1, - 


l\\ 

x4 

> 7L - 
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-^ z - 


^11 

x2 
> 7L 



'^^^ 7ri(i?2)-^7ro(i?3)^7ro(i?)^7ro(i?2) 

^" ^11 .1 



-> 



The Z2-cohomology is then the homology of the hom se- 
quence of this chain complex into Z2: 



H^{Q,dQ;Z2) H\Q,dQ;Z2)^Hi{Q;Z2) 



l\\ 



Poincare duality 



Zo 



















Thus, one finds that H'^{K{Z2, 2); Z2) = Z2, so there are 
only two possible fibrations (up to homotopy type) as above, 
the product K{Z2,2) x ii'(Z2,3) and a "twisted product." 
Fact E.l says that 7r3(T) actually comes by composing the 
Hopf map into the (2-sphere inside the) base -?i'(Z2, 2). This 
is clearly false for the product where that composition would 
factor through tt3{K{Z2, 2)) ^ 1. This shows that the model 
T for U/O has nontrivial k invariant. 

If X is a space and 

F ^ E 

B 

is a fibration, then there is a corresponding fibration of the 
space of maps: 

M{X,F)^M{X,E) 

M{X,B) 

(We now use the usual comma instead of "^" in the no- 
tation for map spaces.) This is our main tool; it enables 
us to tear apart maps into T into maps into ii'(Z2,2) and 
K{Z2, 3) which are easy things to compute, namely cohomol- 
ogy groups. A small wrinkle is that because of the boundary 
conditions we will study relative maps and end up with rela- 
tive cohomology groups, but these are also easily computed. 

So, let us now return to the computation of 7ro(-R2n)- We 
study the fibration: 



The boundary map is always zero for product fibrations, but 
in this case, it is potentially affected by the k invariant (coming 
from the k invariant of T). 

Following the usual procedure, di may be understood as 
the "twist" (mod 2) in a ribbon induced by the "lasso move" 
shown in figure 16 below, which is well known to be 4tt. 




4tt twist 



FIG. 16: By pulling a ribbon over a hedgehog and 'lassoing' it, a 47r 
twist is performed. 

Since Att is an even multiple of 2tt, di is in fact zero and 
7ro(i?) = Z2. The two components can, on the level of framed 
ribbons, be identified with a total twisting of ribbons being an 
even or odd multiple of 2tt. 

TTi (R) may be similarly computed: 



->7ri(i?3) 



-TTliR)- 



-^TTl{R2)- 
l\\ 



H^iQ X /, diQ X /); Z2) H^{Q x /, d{Q x /); Z2) 

^11 < — Poincare duality — >- l\\ 

H^{QxI;Z2) H2{QxI;Z2) 

l\\ l\\ 

1 



q- 



Thus, 7ri(i?) = Z2 ■ Intuitively, 7ri(i?) is represented by 
2-dimensional "bags" or homology classes that a loop of rib- 
bons, defining an element of tti (R), sweeps out in Q x /. 



R3^ R 

R2 

where R is our original R2n (the subscript n > is sup- 
pressed). R2 is the space Maps(Q -^ K{'Z2,2)) with rigid 
boundary conditions (from dQ to S*^ C K{'Z2, 2)) analogous 



Appendix F: Pontryagin-Thom Construction 



Homotopy classes [M'^, S''] may be studied geometrically. 
The idea is to make any map / : Af* -^ S'^ transverse to a 
base point A^ e S''. Then f^^{N) is a (d — A:)-dimensional 
submanifold K of M"^ equipped with a framing of its normal 
bundle obtained by pulling back a fixed normal /c -frame to 
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N. This construction is reversible: Given a submanifold with 
normal framing K'^^^ C M, there is a map / : Af* — > S^ 
which wraps an e-neighborhood rj of K over S'^ by using polar 
coordinates around the north pole N and the framing to send 
each normal D^ C r] degree = 1 over S'' and sends the rest 
of M \ 77 to the south pole in S'^ . 

This discussion can be relativized. A homotopy F : M x 
\0_,1] ^ S, F\mxo ^ f, F\mxi = g yields by inverse image 
K of N a "framed cobordism" from Kq C A/ x to Ki C 

M X 1, that is, a manifold K C M x [0, 1] with a normal 

fc-framing which restricts at (1) to ii^o(i) ^nd its normal k- 
framing. 

The maps to S'' that are produced from a framed submani- 
fold may seem extreme and unrepresentative since most points 
map to the south pole. However, the beauty of the construc- 
tion is that since S'^ \ N is contractible, no real choice (up 
to homotopy) exists for the part of the map which avoids N. 
Thus we have: 

Fact F.l. The space of maps Ai (Af*, S'^) is (at least weakly) 
homotopy equivalent to the space of framed submanifolds 
S^j^d-k ^ TVf^} (provided both are given reasonable topolo- 
gies). 

Fact F.l. Fact F.l holds in a relative setting, e.g. the 

weakly 

{space of homotopies] ~ {space of framed cobordisms). 
For example, merely on the level ofirg, this says that "homo- 
topy classes of maps (Af , S ) are in 7-7 correspondence with 
framed cobordism classes of framed submanifolds K' 

M''." 



d-k 



C 



This basic PTC may be generalized to maps [N'^^M^) 
whenever we can identify a "spine" X C Af*^ so that 
Af \ X is contractible. One then studies f^^{X) (actually 
/^^ (a germ of X) corresponding to the old framing data) in- 
stead of /^^ (iV). For this approach to be practical, X and its 
neighborhood have to be fairly simple. In this extended set- 
ting, one must also keep track of the maps /| : {f^^{X)) -^ 
X since this is no longer unique. 



Appendix G: Multiplication Table for iTi{K2n) 

We have shown 72n — T^i{K2n) = Z x Z2 x 6*2,1 and a 
short exact sequence (SES): 



total parity 

1 >G2n ^Z2"X52„ ^^2 



Thus, we have another SES: 



-> 1 



{a,g){b,h) = (a + g.b + f{g,h),gh) 



(Gl) 



where a e 



^^" g, h e S2n- 9 acts on Z2" ^ by including 



Z2" ^ C Z2" as the "even" vectors and letting g permute the 
factors of Z2". It only remains to determine a twisted cocycle 
/ e H^(S2n;^l"~'^), where the action of 52„ on Z^"-i is 
as above. The cohomology class represented by / is non-zero 
because the exact sequence does not split. As a check, we 
know (line 3.10 in^*") that / should obey: 



g{.f{h, k)) - f(gh, k) + fig, hk) - fig, k) - (G2) 

Claim. fig,h) — if g or h is even and fig,h) ~ 6g(i) 
if g and h are both odd, where eg(i) is the .g(l)''' standard 
generator ofZ'^". 



Appendix H: Ghostly recollection of the braid group 

To promote the Teo-Kane representation to a linear repre- 
sentation of the braid group Bm, we have the freedom to intro- 
duce an overall abelian phase to each Tij. There are several 
different choices related to different physical systems. For our 
purpose here, any choice suffices. The different choices lead 
to different linear images of the braid groups as the braid gen- 
erator matrix has different orders. 

The Ising TQFT can be realized using the Kauffman 
bracket by choosing A = ie^ t^ . The braid group represen- 
tation from a TQFT is not well-defined as a matrix represen- 
tation because in general we do not have a canonical choice 
of the basis vectors. In the Kauffman bracket formulation 
of Ising TQFT, bases of representation spaces can be con- 
structed using linear combinations of Temperley-Lieb (TL) 
diagrams or trivalent graphs. Then the representation of the 
braid groups is given by the Kauffman bracket interpreted as 
a map from Bm to units of TL algebra. 

The Temperley-Lieb-Jones algebra TLmiA) at the above 
chosen A is isomorphic to the Clifford algebra. Therefore the 
7-matrices {ji} can be represented by TL generators {Ui}. In 
terms of the 7-matrices, the Jones representation in Kauffman 
bracket normalization is: 



PAio-t) 



- 4 7i7i+l 



Jones original representation from von Neumann algebra is: 

pi(a,) = eTe-I'''T-+i. 



-^Z 



2n-l 



-^G 



2n 



with a merely set theoretic section s. However, using s (line 
3.5 of^*") we can write the multiplication table out for G'2n, 
hence 7TiiK2n) in terms of a twisted 2-cocycle / on S2n'- 



We will also refer to the 7i i+i representation of the braid 
group 

as the 7-matrix representation. 

The ratio of the two distinct eigenvalues of the braid gen- 
erator is —i, independent of the overall abelian phase. The 
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orders of the braid generator matrices in the three normaliza- 
tions are 16, 4, 8, respectively. In the following, we will focus 
on the Jones representation. 

Let q = A^'^ = i, then the Kauffman bracket is pA {<Ji) =< 
CFi >= A + A^^Ui, where C/i's are the TL generators. For 
convenience, we introduce the Jones generators of the TL al- 
gebras ei = ^ , d ~ \/2, i = 1, • • • ,771—1. These e/s should 
not be confused with the basis element e^'s of M™. The TL 
algebra in terms of e/s is: 



e- e^ 



_ 1 



-i+(i+ 



The Jones representation is p^(cri) — —ApA{(Ji) = - 
q)ei. Let Xi —< af >. Then we have: xf — 1, XiXj — ^j^ 
if \i — j\ > 1, and XiXi+i — — x^+ix^. The —1 in the relation 
XiXi+i = —Xi+iXi is from the Jones-Wenzl projectors p^ = 
0. By a simple calculation, the Jones-Wenzl projector p^ in 
termsof {ei,ei+i}isp3 = l + 2(eiej+i+ei+iej-ei-ei+i), 
and x^Xi+i + Xi+iXi = 2p^. 

Note that {xi = 1 — 2ei,i — I,--- ,to — 1} gener- 
ate TLm{A). Write Xi in terms of 7-matrices, we have 
Xi — a/— T7i7i+i, which identifies TLm{A) with the even 
part C/"„(C) of Clm{C). On the other hand, if we set 



■■Xl,l 



!,■ 



,t 



1, a direct com- 

26, 



, m 
putation shows that v^ — Vi,ViVj + VjVi — 2Sij. Therefore, 
{vi, i = 1, • • • ,m — 1} form the Clifford algebra Clm-ii'C). 
Of course, the two different realizations TL„i{A) as Clifford 
algebras are just the well known isomorphism C/J'„(C) = 

Clm-l'C). 

To identify the image of the Jones representation, we use 
the exact sequence 1 — > PBm — > Bm -^ Sm -^ 1- The im- 
age of the pure braid group PB^ is generated by {x^}. As 
an abstract group, the group can be presented with generators 

Xi,i = !■■■ ,Xm-i and relations xf = 1, XiXj — XjXi if 



Xi+iXi. This group, called the 
^'^ and denoted as E}^_i, 



is of 



\i- j\ > 1, andXtXi+i - 
nearly-extra-special 2-group in 

order 2™™!. Consequently, the image Gm of the Jones rep- 
resentation of Bm fits into the exact sequence 1 -^ Em-i ^ 
Gm — > Sm — > 1. If we projectivize this sequence, we have 



1 



I/J, 



Gn 



bn 



1 . (This sequence splits if and 



only if m is odd.) In this sense, the ribbon permutation group 
is a ghostly recollection of the braid group. 

The center Z{E^_i) of £'„_! is Z2 = {±1} if "t- is even, 
and Z2 X Z2 = {±l,±xiX3 • • ■ Xm-i} if rn is odd. The 
representation of Ejj^_^ is faithful, but when m is even, none 
of the irreducible sector is faithful because the central element 
zLxiXs ■ ■ ■ Xm-i acts by ±1, too. The action of Xi on the two 
irreducible sectors are the same for i ~ 1,2,--- ,m — 2 and 
differs only on Xm-i, which is ±1. When projectivized, note 
that the element a = (11 • • • 1) e Z™, as an element in Z™"\ 
is invariant under the action of Sm, hence the projective image 
of each irreducible sector is isomorphic to Z™^^. 

As comparison, we mention the related result in^^. The 
Jones representation there is in the 7-matrix normalization. 
The pure braid generators are of order 4. Hence the pure braid 
group image is a variant -E^L^ of E}^_i. The generators of 
the two groups are related by the change of Xi to —y^—lxi. 
A presentation of E^^_^ with generators Xi,i = 1 ■ ■ ■ , Xm-i 
has relations xf = —1, XiXj = XjXi if \i — j\ > 1, and 
XiXi+i = —Xi+iXi. The center Z{E^_-^) is more compli- 
cated: Z2 if m is odd, and Z2 x Z2 if m = 4fc + 1 and Z4 if 
TO = 4fc + 3. When to. is even, the image of each irreducible 
sector is identified as the lifting of Gm to Spin{m)^^ . It fol- 
lows from the above discussion, the linear image of the Jones 
representation in 7 matrix-normalization fits into the exact se- 
quence 1 ^- Z4 X El^_j^ — > Zg X Gm ^ Z2 X S*™ — > 1, 

where Gm is the image of Bm above. 
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